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CALCULATION  OF  VISCOUS  EFFECTS  ON  SHIP  WAVE  RESISTANCE 
USING  AXISYMMETRIC  BOUNDARY  LAYER  APPROACHES 


INTRODUCTION 

An  historic  paper  by  Micheli  [1]  developed  a  direct  calculation  procedure  based  on  thin-ship 
theory.  Since  then,  calculation  of  the  resistance  due  to  the  waves  generated  by  a  surface  ship  advancing 
at  constant  forward  speed  has  been  an  area  of  active  interest.  Several  international  seminars  and 
workshops  on  wave  resistance  have  been  held  in  recent  years.  The  accurate  calculation  of  the  wave 
resistance  is  important  for  commercial  as  well  as  Navy  ships  since  it  is  one  of  the  major  components  of 
overall  drag.  In  addition,  since  the  wave  resistance  represents  a  particular  weighted  integral  of  the 
amplitude  spectrum  of  the  waves,  its  accuracy  also  gives  an  indirect  indication  of  the  overall  accuracy  of 
the  calculated  shape  of  the  wave  spectrum  in  wave  number  space  and  the  associated  wave  pattern  on 
the  physical  free  surface. 

Most  of  the  efforts  aimed  at  improving  the  agreement  of  the  calculated  and  measured  results  have 
centered  along  two  distinct  lines.  In  one  approach,  potential  flow  is  assumed;  efforts  are  made  to  better 
satisfy  the  kinematic  condition  on  the  hull  surface  and  the  nonlinear  dynamic  condition  on  the  un¬ 
known  free  surface.  Due  to  the  complexity  of  free  surface  flows,  these  efforts  often  lead  to  calculation 
procedures  which  are  orders  of  magnitude  more  time-consuming  than  the  original  Micheli  approach.  In 
the  other  approach,  the  assumption  of  potential  flow  is  removed,  and  attempts  are  made  to  account  for 
the  effect  of  fluid  viscosity  on  the  flow  near  the  ship.  The  usual  procedure  is  to  take  the  viscous  flow 
to  effectively  change  the  shape  of  the  ship  in  the  stern  and  near-wake  regions.  The  proposed  calcula¬ 
tion  procedures  have  ranged  from  the  use  of  empirical  correction  factors  or  the  assumption  of  plausible 
mathematical  shapes  to  the  use  of  complex  and  time-consuming  higher  order  boundary  layer  theory. 
The  effect  of  the  resultant  change  in  hull  shape  on  the  calculated  wave  resistance  is  often  assessed  by 
using  the  original  thin-ship  theory. 

This  report  investigates  the  effect  of  fluid  viscosity  for  a  scope  which  is  wider  than  that  of  previ¬ 
ous  studies.  Two  direct,  simple  models  are  used  to  model  the  potential  flow:  the  previously  mentioned 
Micheli  thin-ship  theory,  and  zeroth  order  slender-ship  theory,  as  it  is  termed  by  Noblesse  [2].  The 
viscous  flow  is  modeled  by  two  relatively  simple,  integral  boundary  layer  methods,  both  of  which 
approximate  the  actual  hull  by  axisymmetric  bodies.  One  method  tends  to  overestimate  the  viscous 
effects  while  the  other  tends  to  underestimate  these  effects.  The  various  methods  are  used  to  calculate 
wave  resistance  for  a  series  of  speeds  for  five  ship  hulls,  three  of  which  have  extensive  experimental 
data.  The  fourth  hull  illustrates  a  case  where  the  present  simple  methods  basically  break  down;  the 
fifth  hull  illustrates  the  effects  of  viscosity  for  an  unsymmetrical  ship  moving  in  both  forward  and  back¬ 
ward  directions.  The  entire  calculation  procedure  is  simple  enough  to  run  with  ease  on  a  minicom¬ 
puter. 

The  report  starts  with  a  statement  of  the  formulation  of  the  wave  resistance  problem  and  a  brief 
review  of  the  various  calculation  methods  used  in  the  potential  and  viscous  flow  areas.  The  review  is 
not  comprehensive  but  indicates  the  principal  approaches  in  each  area.  A  description  is  next  given  of 
the  relatively  simple  potential  flow  and  axisymmetric  boundary  layer  methods  used  in  the  present  study. 
Their  limitations  and  likely  accuracy  are  discussed.  The  major  geometrical  features  of  the  five  ship 
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hulls  considered  in  the  present  study  are  given.  The  calculated  boundary  layer  parameters,  as  well  as 
the  calculated  and  measured  values  of  wave  resistance,  are  presented  in  tabular  and  graphic  forms.  The 
calculated  values  of  wave  resistance  with  viscous  correction  effects  are  compared  with  corresponding 
potential  now  results  and  with  measured  values.  The  report  concludes  with  a  summary  of  the  principal 
results. 


FORMULATION  OF  THE  WAVE  RESISTANCE  PROBLEM 

The  problem  is  usually  formulated  for  a  coordinate  system  fixed  to  the  ship.  Figure  1  shows  the 
coordinate  system  used  in  the  present  study,  with  origin  at  the  forward  perpendicular,  x  positive  from 
stern  to  bow,  y  positive  to  port,  and  z  positive  upwards.  The  figure  also  defines  terms  for  the  ship 
configuration  which  will  be  used  throughout  this  report.  In  this  coordinate  system,  the  flow  pattern  is 
steady  with  uniform  flow  velocity  U  in  the  negative  x -direction  as  x  -+  +  oo.  The  assumptions  are 
usually  made  that  the  fluid  is  inviscid  and  incompressible,  the  flow  is  irrotational,  and  surface  tension 
effects  may  be  neglected.  Refinements  to  account  for  viscous  effects,  which  are  important  in  the  stern 
and  near  wake,  are  discussed  later.  Surface  tension  effects  are  important  only  for  waves  of  small 
wavelength,  such  as  capillary  waves. 


Fig.  1  —  Ship  configuration  and  coordinate  system 


Under  the  above  assumptions,  the  flow  field  is  described  by  the  velocity  potential  <Mxj>,z)  which 
satisfies  Laplace’s  equation 

V2  d>  =  0.  (1) 

On  the  free  surface  z  =  the  dynamic  and  kinematic  boundary  conditions  are  respectively  given 

by 

g  C  +  y(v<i>)2=  y  U1  (2) 

and 


■Jt  (z  -  0  =  V  <&  •  V  (z  -  0  =  <t>Ax  +  -$z  =  0  (3) 

where  g  is  the  gravity  constant  and  t  is  time.  The  combination  of  the  dynamic  and  kinematic  condi¬ 
tions  gives 


g<Fz  +  V  $>  •  V 


y(V<F)2 


=  0 


on  2  =  £(x,y). 


(4) 
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The  kinematic  boundary  condition  on  the  wetted  ship  hull  SH  is  that  no  fluid  goes  through  the  ship  hull 

n  •  V4>  =  0  on  5^  (5) 

where  n  is  the  inward  normal  to  the  ship  hull.  Finally,  the  radiation  condition  states  that  waves  are  not 
radiated  upstream. 

Once  <J>  is  obtained,  the  wave  resistance  Dw  may  be  obtained  by  direct  integration  of  the  pressure 
p  acting  on  the  hull  surface 

Dw  =  /  Js  pnx  ds  (6) 

where  nx  is  the  x-component  of  the  inward  normal  to  the  hull,  and  p  is  related  to  by  Bernoulli’s 
equation 

P  =  -  -^[(V*)2-  U2]  -pgz  (7) 

where  p  is  the  fluid  density. 

In  the  above  form,  the  calculation  of  Dw  requires  the  evaluation  of  gradients  of  <t>,  which  tends 
to  accentuate  inaccuracies  in  O.  An  alternate,  more  commonly  used,  approach  is  to  relate  the  wave 
resistance  to  the  energy  flux  of  the  generated  waves  across  a  vertical  plane,  x  =  constant,  downstream 
of  the  ship.  Following  the  derivation  given,  for  example  in  Refs.  3  or  4,  the  following  expression  for 
Dw  in  terms  of  the  amplitude  spectrum  A2 (9)  is  obtained 

Dw  =  —  TT  p  U 2  J_v/2  A2 (9)  cos3  9  d9  (8) 

where  9  is  the  angle  of  direction  of  propagation  with  respect  to  the  x-axis  of  the  component  wave.  For 
the  usual  case  where  $  is  represented  by  a  distribution  of  sources  over  the  hull  surface,  it  is  further 
shown  in  Ref.  4  that  A2 (9)  is  given  by  an  integral  of  the  source  densities  over  the  hull  surface  weighted 
by  a  wave  elevation  function. 

Equation  (8)  shows  that  the  wave  resistance  represents  an  integrated  average  of  the  amplitude 
spectrum  weighted  by  the  factor  cos3  9.  This  means  that  the  wave  resistance  is  dominated  by  the 
transverse  waves,  for  which  |#|  <  35°.  Using  the  relations  given  in  Ref.  4,  the  amplitude  spectrum  can 
be  alternatively  expressed  as  a  single  function  of  the  wavenumbers  in  the  xory  directions,  respectively 
kx  and  ky.  These  relations  in  terms  of  kx  or  ky  would  be  useful  in  obtaining  wave  resistance  from  a 
longitudinal  or  lateral  wave  cut. 

PREVIOUS  STUDIES 

Potential  Flow 

The  Proceedings  of  the  Workshop  on  Ship  Wave-Resistance  Computations  [5]  gives  numerous 
examples  of  recent  approaches  for  solving  the  above  problem.  Relatively  few  studies  attempt  to 
directly  solve  the  problem,  which  includes  the  application  of  the  nonlinear  free  surface  condition  (Eq. 
(4))  on  the  unknown  free  surface  z  —  £.  Chan  and  Chan  [6,7]  obtain  the  steady  state  solution  by 
marching  forward  in  time  and  using  a  finite  difference  grid  to  model  the  flow  field.  Their  approach  is 
unsuitable  for  the  calculation  of  the  wave  resistance  since  the  solution  must  be  solved  for  a  number  of 
time  steps  until  the  steady  state  solution  is  reached.  Von  Kerczek  and  Salvesen  [8-10]  also  use  a  Finite 
difference  approach  to  solve  the  flow  around  prescribed  isolated  and  distributed  singularities.  However, 
their  work  is  restricted  to  two-dimensional  flow.  At  a  somewhat  lower  level  of  complexity,  coordinate 
transformation  techniques  may  be  used  to  transform  the  physical  coordinate  system  into  a  new  coordi¬ 
nate  system  where  the  free  surface  is  known.  The  actual  nonlinear  problem  is  converted  to  a  series  of 
linearized  problems  on  the  successively  iterated  free  surface.  Two  methods  which  have  been  used  are 
Guilloton’s  method  [11,12]  where  lines  of  constant  pressure  become  constant  coordinate  lines,  and  a 
Lagrangian  method  [13]  where  streamlines  become  constant  coordinate  lines. 
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The  majority  of  studies  consider  only  the  linearized  form  of  Eq.  (4)  and  apply  this  condition  oil 
the  undisturbed  free  surface  z  =  0  instead  of  on  z  =  £.  Since  the  disturbances  to  the  free  surface  are 
considered  to  be  small,  it  is  convenient  to  express  the  total  potential  $  as  the  sum  of  a  base  flow  poten¬ 
tial  <f>B  and  a  perturbation  potential  0;  here 

<l>  =  <f>B  +  0.  (9) 

The  base  flow  is  usually  taken  to  be  uniform  flow,  for  which 

0B  =  -Ux,  (10a)  4 

or  the  flow  around  the  ship  hull  where  the  free  surface  is  taken  to  be  a  horizontal  plane  of  symmetry, 
for  which 

5  Si  J 

0b  =  4>d  ■  (10b) 

The  latter  case  is  often  referred  to  as  double  model  flow  and  may  be  viewed  as  the  limiting  case  of 
velocities  U  which  are  slow  enough  that  the  free  surface  is  undisturbed.  If  the  base  flow  potential  is 
given  by  Eq.  (10a),  the  linearized  form  of  the  free  surface  condition,  Eq.  (4),  is 

g<f>z  +  U2^  =  0  (ID 

and  the  hull  boundary  condition,  Eq.  (5),  becomes 

n  ■  Vd>  =  Unx.  (12) 

Two  computer  methods  developed  by  Chang  [14]  and  Dawson  [15]  are  representative  of  the 
approaches  used  to  solve  the  above  (linearized)  problem.  Chang  divides  the  ship  into  a  series  of  quad¬ 
rilaterals,  on  each  of  which  is  placed  a  Kelvin  (or  Havelock)  source  of  uniform  but  unknown  density  a. 
The  Green  function  for  the  Kelvin  source  G  [16,17]  may  be  written  as 

G  =  —  -  —  +  G2.  (13) 

r  r  l 

The  first  term  is  simply  the  potential  for  a  Rankine  (infinite  fluid)  source  placed  on  the  ship  hull 
(x0,y0,  -  z0).  The  second  term  is  the  potential  for  a  Rankine  sink  placed  at  (x0,y05  +  zo)-  The  third 
term  is  a  double  integral  which  gives  the  wave  disturbance  due  to  an  oscillating  source  placed  below  a 
linearized  free  surface.  Since  the  Kelvin  source  is  specifically  constructed  to  satisfy  the  linearized  free 
surface  condition  (Eq.  (11))  and  the  radiation  condition,  the  hull  boundary  condition  (Eq.  (12))  results 
in  a  series  of  algebraic  equations  for  the  unknown  source  strengths  cr .  The  principal  difficulty  in  this 
approach  is  the  careful  evaluation  of  the  matrix  of  velocity  influence  coefficients  due  to  the  complex 
Havelock  sources.  Also,  matrix  inversion  techniques  must  be  used  to  solve  the  resulting  set  of  alge¬ 
braic  equations.  The  type  of  problem  solved  by  the  Chang  approach  is  commonly  called  the 
Neumann-Kelvin  problem. 

Dawson  uses  a  different  approach  which  avoids  the  above  difficulty  involving  the  use  of  Havelock 
sources.  He  linearizes  the  free  surface  condition  (Eq.  (4))  by  using  the  double  model  flow  as  the 
known  base  flow,  resulting  in 

0 0/0 Dll  +  (0O  0/)/  +  S  4>z  =  0  ^) 

where  /  denotes  differentiation  along  a  streamline  of  the  flow.  He  divides  the  actual  ship  hull,  its 
reflection  above  the  undisturbed  free  surface,  and  the  undisturbed  free  surface  adjacent  to  the  hull  into 
a  series  of  quadrilaterals.  He  places  a  Rankine  (infinite  fluid)  source  of  uniform  but  unknown  density 
on  each  quadrilateral.  The  kinematic  hull  condition  (Eq.  (5))  is  used  for  the  quadrilaterals  on  the  hull 
and  its  reflection.  The  free  surface  condition  (Eq.  (14))  expressed  in  finite  difference  form,  is  used  for 
the  quadrilaterals  on  the  free  surface  near  the  hull.  A  four-point  difference  operator  is  found  to 
approximately  satisfy  the  radiation  condition  in  that  upstream  waves  are  suppressed.  The  resulting  set 
of  algebraic  equations  can  then  be  solved  for  the  strength  cr  of  the  source  on  each  quadrilateral.  One  of 
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the  principal  difficulties  of  this  approach  is  that  panels  must  be  used  to  model  not  only  the  actual  hull 
but  also  its  reflection  and  part  of  the  free  surface.  Also,  great  care  must  be  exercised  in  choosing  the 
proper  extent  of  the  free  surface  and  the  manner  of  application  of  the  free  surface  condition. 


The  above  shows  that  the  solution  of  even  the  linearized  problem  is  a  rather  formidable  task. 
The  solution  can  be  simplified  considerably  if  the  beam  B  and  draft  H  are  small  compared  to  the  ship 
length  and  the  dominant  length  U/Vg  of  the  generated  waves  so  that  the  disburbance  potential  <f>  can 
be  considered  negligibly  small  compared  to  the  base  flow  potential  cf>B.  In  this  case,  termed  zeroth 
order  slender-ship  theory  by  Noblesse  [2],  the  source  strength  a  can  be  directly  obtained  from  the 
geometry  of  the  hull 


a  (x,y,z) 


—  U 

+df/dx 

47 T 

Vl  +  (df/dx)2  +  (df/dz)2 

(15) 


where  /6c,z)  is  the  local  half-width  of  the  ship.  The  resulting  direct  expression  for  the  perturbation 
potential  <f>  is  given  by  two  integrals  which  are  entirely  dependent  on  ship  geometry: 

0  =  Js  a-  G  ds  +  F2  cr  nxry(T  G  dl ,  (16) 


where  Fn  =  U/VgL  is  the  Froude  number, 

C  is  the  waterline  contour  of  the  ship,  and 
t y  is  the  ^-component  of  the  tangent  to  the  waterline  contour. 

Noblesse  shows  that  a  series  of  higher  order  slender-ship  solutions  may  be  obtained  by  an  iterative  pro¬ 
cedure  whereby  the  nth  order  solution  is  obtained  as  a  perturbation  to  the  uniform  base  flow  plus  the 
( n  —  l)st  order  solution.  Thus,  the  above  zeroth  order  solution  is  the  perturbation  to  the  base  flow 
only,  neglecting  any  flow  caused  by  motion  of  the  ship.  On  the  other  hand,  the  first  order  solution 
would  be  the  perturbation  to  the  base  flow  plus  the  above  zeroth  order  solution.  In  the  remainder  of 
this  report,  the  zeroth  order  solution  given  by  Eqs.  (15)  and  (16)  will  simply  be  referred  to  as  slender- 
ship  theory. 

A  further  simplification  occurs  if  the  beam  of  the  ship  is  also  assumed  to  be  small  compared  to 
draft,  as  well  as  the  previously  mentioned  ship  and  wave  lengths,  resulting  in  the  Michell  thin-ship 
theory.  In  this  case,  the  surface  integral  in  Eq.  (16)  is  calculated  over  the  centerplane  50  instead  of  the 
actual  hull  surface  Sh>  and  the  contour  integral  is  negligibly  small.  Also,  the  terms  involving  (df/dx)2 
and  (df/dz)2  in  Eq.  (12)  are  neglected.  Equations  (12)  and  (13)  then  respectively  become 


(x,z)  = 

-2  U  df 

(17) 

47 T  dx 

c r  G  ds 

(18) 

where  the  factor  2  in  Eq.  (17)  indicates  that  both  sides  of  the  ship  hull,  assumed  to  have  port  and  star¬ 
board  symmetry,  are  collapsed  onto  the  centerplane. 

Stern  Boundary  Layer  Flow 

The  viscous  boundary  layer  in  the  stern  region  is  usually  of  sufficient  thickness  to  significantly 
influence  the  flow  around  the  ship.  An  indication  of  the  importance  of  the  viscous  flow  is  that  several 
characteristic  differences  exist  between  the  wave  resistance  values  calculated  using  methods  which 
neglect  the  effects  of  viscosity  and  measured  results.  The  largest  difference  is  that  measured  results 
show  reduced  amplitudes  of  variation  with  the  Froude  number  compared  to  theoretical  values.  A 
second  difference  is  that  the  location  of  the  humps  and  hollows  of  the  wave  resistance  curve  for  mea¬ 
sured  values  tend  to  be  shifted  to  higher  Froude  numbers. 
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Qualitatively,  the  boundary  layer  may  be  viewed  as  essentially  thickening  the  width  of  the  ship, 
leading  to  smaller  magnitudes  of  the  slope  S//9x  in  the  stern  region.  This  in  turn  means  a  reduction 
in  the  strength  and  wave-making  ability  of  the  sinks  located  in  the  stern  region.  Thus,  the  earliest 
efforts  at  accounting  for  the  effect  of  viscosity  centered  on  finding  simple  or  empirical  expressions  for 
this  reduction.  Havelock  [18]  found  that  a  constant  reducing  factor  of  f3  =  0.6  for  the  sink  strengths 
over  the  rear  half  of  the  ship  significantly  reduced  the  amplitude  of  oscillation  of  the  wave  resistance 
curve.  Wigley  [19]  deduced  from  experiments  the  following  reducing  factor  1-/3^  for  the  wave  resis¬ 
tance  of  the  after-body  of  the  ship: 

1  -  /3j  =  1  -  e-°m/Fn5.  (19) 

This  shows  that  the  wave  resistance  reducing  factor  1  —  /3^  is  largest  at  small  Fn  and  approaches  0  (i.e., 
no  reduction)  as  Fn  increases. 

Later  efforts  centered  on  the  effect  of  different  assumptions  of  stern  and  near  wake  shapes  on  the 
wave  resistance  curve.  Havelock  [20]  and  Milgram  [21]  make  parametric  studies  of  the  effect  of  dif¬ 
ferent  shape  modifications  in  the  stern— near  wake  region  on  the  wave  resistance  curve.  They  both 
assume  that  the  boundary  layer  is  negligibly  thin  up  to  a  certain  ’’separation'*  point  on  the  aft  body. 
Downstream  of  this. point,  a  mathematical  form  (with  thickness  and  slope  equal  to  that  of  the  actual 
hull  at  the  separation  point)  is  assumed.  Milgram  assumes  a  form  which  has  zero  slope  from  the  aft, 
perpendicular  to  infinity.  This  leaves  an  open  wake  of  constant  width  aft  of  the  ship.  Havelock 
assumes  a  quartic  (or  cubic)  curve  which  has  zero  slope  and  zero  (or  nearly  zero)  thickness  at  various 
selected  distances  behind  the  ship.  This  gives  rise  to  a  wake  of  zero  (or  small)  thickness  beyond  this 
point.  In  both  studies,  the  amplitude  of  oscillation  of  the  wave  resistance  curve  tends  to  decrease  with 
forward  positioning  of  the  separation  point. 

The  above  studies  are  largely  concerned  with  giving  overall  trends  of  the  possible  effect  of  viscous 
flow  on  wave  resistance.  As  such,  they  cannot  be  expected  to  yield  accurate  results  for  a  particular 
ship.  More  recent  studies  have  largely  used  either  low  order  boundary  layer  theory  or  experimental 
measurements  to  obtain  the  displacement  thickness  8*.  This  is  then  added  to  the  stern  and  near  wake 
region  of  the  ship,  resulting  in  an  effectively  thicker  and  longer  ship.  The  earliest  attempts  along  these 
lines  are  those  by  Wigley  [22]  and  Wu  [23]  who  use  two-dimensional  boundary  layer  theory  to  calculate 
8*  up  to  the  aft  perpendicular.  For  the  wake  region,  Wigley  uses  an  approach  similar  to  the  previously 
mentioned  Havelock  [20]  approach  to  obtain  a  closed  (or  nearly  closed)  wake  at  a  certain  distance 
behind  the  ship.  He  shows  that  the  calculated  wave  resistance  is  sensitive  to  the  wake  distance  but 
gives  no  guidelines  for  its  choice.  Wu  mentions  various  ways  of  treating  the  wake  region  but  essentially 
leaves  it  out  in  his  numerical  examples. 

Himeno  [24],  Kinoshita  [25],  and  Larsson  and  Chang  [26]  make  rather  comprehensive  studies  of 
the  effect  of  different  boundary  layer  and  wake  approximations  on  the  calculated  wave  resistance  for 
particular  ship  forms.  Himeno  uses  three-dimensional  integral  boundary  layer  theory  to  calculate  the 
tangent  and  cross  flows  along  the  streamlines.  He  uses  Michell  thin-ship  theory  to  calculate  the  wave 
resistance,  with  and  without  viscous  effects,  for  the  Series  60,  CB  =  0.70  form.  He  finds  that  the  wave 
resistance  curve,  for  the  case  where  only  the  displacement  thickness  for  tangential  flow  is  included, 
gives  an  appreciable  attenuation  of  the  hollows  and  bumps  at  low  Froude  numbers.  However,  the  wave 
resistance  curve  for  the  more  complete  case  where  displacement  thicknesses  for  both  tangential  and 
cross  flows  are  included,  has  pronounced  hollows  and  humps  and  is  in  general  higher  than  the 
corresponding  curve  with  no  viscous  corrections.  This,  of  course,  is  contrary  to  the  trend  for  most 
measured  cases. 

Kinoshita  [25]  gives  an  estimate  of  the  various  effects  of  viscous  flow  on  wave  resistance.  His 
main  conclusion  is  that  the  effect  on  wave  propagation  and  the  effect  due  to  the  rotational  flow  of  the 
wake  are  small  compared  to  the  deformation  of  the  potential  flow,  as  modeled  by  the  displacement 
thickness.  He  uses  measured  boundary  layer  velocity  distributions  to  obtain  the  displacement  thickness 
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along  the  hull  and  in  the  wake.  His  results  for  a  thin  strut  show  that  the  inclusion  of  the  measured  dis¬ 
placement  thicknesses  moves  the  Michell  wave  resistance  curve  into  substantially  good  agreement  with 
measured  values.  His  results  also  show  the  interesting  trend  that  there  is  relatively  little  difference  in 
the  wave-resistance  values  for  the  case  where  8*  is  assumed  to  be  constant  with  depth  and  the  more 
realistic  case  where  8  *  varies  with  depth. 

Larsson  and  Chang  [26]  use  higher-order  boundary  layer  theory  to  calculate  the  displacement 
thickness  along  the  hull  and  in  the  near  wake  for  the  Series  60,  CB  =  0.70  ship  form  used  by  Himeno 
[24],  They  make  wave-resistance  calculations  using  Michell  thin-ship  theory  as  well  as  the  method  by 
Chang  [14],  using  the  more  complex  Neumann-Kelvin  approach  given  by  Eqs.  (9)-(13).  For  the  thin- 
ship  case,  they  make  calculations  for  the  wake  terminating  at  the  aft  perpendicular  (leading  to  a  rela¬ 
tively  large  open  wake)  as  well  as  using  a  calculated  wake;  this  results  in  a  narrower  but  still  open  dis¬ 
placement  body  aft  of  the  ship.  Even  though  no  experimental  data  are  shown,  several  important  points 
may  be  deduced  from  their  calculated  results.  First,  there  is  relatively  good  agreement  between  the  cal¬ 
culated  results,  neglecting  viscous  effects,  using  the  thin-ship  and  Neumann-Kelvin  approaches. 
Secondly,  the  effect  of  the  viscous  correction  is  to  lower  the  wave  resistance  in  all  cases.  Thirdly,  the 
reduction  in  wave  resistance  is  greater  for  the  narrower  wake,  but  the  difference  in  the  reduction 
between  the  two  wakes  varies  with  Froude  number.  They  show  that  thin-ship  theory,  when  applied  to 
an  open  wake,  does  not  properly  account  for  the  contribution  of  the  transverse  wave  system  to  the 
wave  resistance.  Thus,  when  the  drag  due  to  the  transverse  waves  is  dominant,  there  is  a  large  differ¬ 
ence  in  the  calculated  reduction  between  the  two  wakes.  On  the  other  hand,  the  difference  is  small 
when  the  divergent  waves  are  dominant.  It  is  of  interest  to  note  that  in  spite  of  their  observation  con¬ 
cerning  the  open  wake,  their  viscous  flow  calculation  procedure  still  leaves  an  open  wake  aft  of  the 
ship. 

PRESENT  CALCULATION  PROCEDURES 
Potential  Flow  Calculation 

Largely  due  to  the  simplicity  of  the  direct  solution  methods,  it  was  decided  to  use  slender-ship 
and  Michell  thin-ship  theories  to  make  the  wave  resistance  calculations.  These  theories  are  respectively 
given  by  Eqs.  (15)-(16)  and  Eqs.  (17)-(18).  Use  of  these  simple  theories  allows  calculation  of  wave 
resistance,  with  and  without  viscous  correction,  for  a  wider  range  of  ship  hull  cases  than  in  previous 
studies.  Also,  the  preceding  discussion  of  previous  studies  shows  that  in  many  cases  the  viscous  flow 
correction  is  as  important  as  higher  order  potential  flow  theory  in  bringing  the  calculated  results  into 
better  agreement  with  measured  values. 

The  thin-ship  calculations  are  made  using  a  computer  program  developed  by  Chapman  [27],  based 
on  the  method  outlined  in  Ref.  4.  Briefly,  the  ship  centerplane  is  divided  into  a  series  of  quadrilaterals 
or  triangles,  on  each  of  which  is  placed  a  source  of  constant  strength  cr,  given  by  Eq.  (17).  The  pro¬ 
gram  then  computes  an  amplitude  spectrum  A 2  for  the  far-field  waves  as  a  function  of  the  longitudinal 
wave  number  kx  by  integrating  the  product  of  cr  and  essentially  the  far-field  wave  part  of  the  Kelvin 
source  G  given  in  Eq.  (13).  The  wave  resistance  is  then  calculated  by  using  an  equation  similar  to  Eq. 
(8),  with  the  exception  that  the  independent  variable  is  now  kx  instead  of  0.  Reference  4  shows  that 
the  amplitude  spectrum  can  be  expressed  as  a  single  function  in  terms  of  kXJ  0,  or  ky,  the  lateral  wave 
number.  It  may  also  be  noted  that  the  wave  elevation  on  the  free  surface,  £(x,y),  may  be  conveniently 
obtained  from  the  amplitude  spectrum  by  using  Fourier  transform  methods. 

Two  extensions  were  added  to  the  program  to  upgrade  it  to  slender-ship  capability.  First,  the  cal¬ 
culation  of  the  line  integral  in  Eq.  (16)  was  added  to  the  program  in  the  form  of  a  subroutine. 
Secondly,  a  geometry  subroutine  [28]  was  added  to  calculate  average  values  of  nx  in  Eq.  (15)  over  a 
ship  panel.  An  additional  consideration  is  that  the  ship  panels  are  now  placed  over  the  actual  hull  sur¬ 
face  SHy  instead  of  the  centerplane  S0,  since  the  surface  integral  in  Eq.  (16)  is  now  over  SH. 
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The  modifications  have  been  made  in  such  a  manner  that  the  use  of  thin-ship  or  slender-ship 
theory  is  at  the  option  of  the  user.  The  program  also  provides  the  option  of  modifying  the  hull  surface 
to  account  for  the  displacement  thickness  which  is  calculated  by  the  boundary  layer  methods  described 
below. 

Boundary  Layer  Calculation 

The  results  of  the  previous  studies  indicate  two  principal  points.  First,  the  wave  resistance  is  sen¬ 
sitive  to  the  shape  of  the  wake  aft  of  the  ship.  The  discussion  in  Ref.  26  indicates  that  an  open  wake 
will  tend  to  overestimate  the  wave  resistance,  with  the  error  depending  on  the  relative  importance  of 
the  transverse  waves.  However,  the  results  of  Ref.  22  indicate  that  closing  the  wake  arbitrarily  may 
lead  to  increases  in  wave  resistance.  Except  for  those  cases  where  the  wake  is  mathematically  shaped, 
the  wake  usually,  has  been  open.  Secondly,  the  results  of  Refs.  24  and  25  indicate  that  more  careful 
modeling  or  calculation  of  the  boundary  layer  effects  need  not  lead  to  wave  resistance  values  which  are 
in  better  agreement  with  experimental  results. 

In  view  of  the  above,  and  the  simplicity  of  the  methods  chosen  to  make  the  wave  resistance  cal¬ 
culations,  it  was  decided  to  use  a  relatively  simple,  axisymmetric  momentum-integral  boundary  layer 
method  to  model  the  overall  characteristics  of  the  viscous  flow  on  the  ship  hull.  Following  the  sugges¬ 
tion  by  Granville  [29],  the  actual  hull  is  approximated  by  an  axisymmetric  body,  as  described  below. 
This  axisymmetric  approach  does  not  appear  to  have  been  used  in  previous  studies  which  investigated 
the  effect  of  viscosity  on  wave  resistance.  These  studies  have  largely  used  two-  and  three-dimensional 
theories.  Because  of  the  complexity  and  uncertainty  surrounding  the  near-wake  flow,  the  boundary 
layer  calculations  are  terminated  at  the  aft  perpendicular.  This  leaves  an  open  wake  which  may  overes¬ 
timate  the  wave  resistance  (i.e.,  underestimate  the  reduction  due  to  viscous  flow)  for  certain  Froude 
number  cases. 

Equivalent  Axisymmetric  Bodies 

In  the  present  approach,  the  three-dimensional  hull  is  replaced  by  an  equivalent  axisymmetric 
body.  Two  possibilities  exist:  respectively  equal  in  volume,  and  wetted  area  to  a  double  model  of  the 
hull,  which  is  defined  as  the  underwater  portion  of  the  ship  and  its  reflection  about  the  free  surface, 
taken  to  be  horizontal.  At  a  given  longitudinal  station  x,  the  radii  of  the  equal  volume  and  equal  wet¬ 
ted  area  bodies,  respectively  Rv  and  Rw,  are  given  by  Ref.  29  as 


(20) 


and 


(21) 


where  A  ix)  is  the  cross-sectional  area  of  the  underwater  portion  of  the  ship,  and 
P(x)  is  the  perimeter  of  the  underwater  portion  of  the  ship. 

Following  Ref.  29,  the  pressure  field  is  calculated  by  using  the  Rv  body  since  it  models  the  volume 
change  along  the  actual  hull.  The  boundary  layer  calculations  are  carried  out  using  both  bodies.  It  is 
pointed  out  later  that  the  use  of  the  two  bodies  tends  to  bracket  the  effect  of  streamline  convergence 
obtained  by  more  complex  boundary  layer  methods. 

Momentum  Integral  Equation 

The  momentum  integral  equation  for  boundary  layer  flow  along  an  axisymmetric  body  may  be 
conveniently  derived  from  the  more  general  three-dimensional  formulation  for  flow  along  the  ship  hull. 
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In  the  more  general  approach,  typically  a  series  of  streamlines  are  calculated  along  the  hull  surface. 
Then  two  momentum  integral  equations  are  integrated  for  tangential  flow  and  cross  flow  along  each 
streamline.  If  one  assumes  all  cross  flow  terms,  as  well  as  derivatives  normal  to  the  streamlines,  as 
being  of  higher  order  and  considers  the  equation  for  tangential  flow  for  a  single  representative  stream¬ 
line,  the  axisymmetric  momentum  integral  equation  [30]  is  obtained 


dO  Cf  ,  |  g  dUe  Q  dR 
ds  2  Ue  ds  R  ds' 


(22) 


Here, 


=r 

J  0 


1  - 


Ug 


——  dX  is  the  momentum  thickness, 

Up 


8  is  the  boundary  layer  thickness  where  -j—  —  0.995, 

V  p 


u  is  the  velocity  in  the  boundary  layer, 

Ue  is  the  velocity  at  the  edge  of  the  boundary  layer, 
s  is  the  arc  length  measured  on  the  body,  starting  from  the  bow, 

Cf  =  Tw/-j-p  Ug  is  the  skin-friction  coefficient, 


tw  is  the  shear  stress  on  the  body  surface, 
H  =  8  */9  is  the  shape  parameter, 


-f 

Jo 


1 


u 

Ug 


dX  is  the  displacement  thickness,  and 


R  is  the  radius  of  the  body. 


(23a) 


(23b) 

(23c) 

(23d) 


Figure  2  shows  the  coordinate  system  used  in  the  boundary  layer  calculations,  where  X  is  positive  from 
bow  to  stern,  contrary  to  the  sense  for  x  for  the  ship  hull  in  Fig.  1.  (This  is  because  the  boundary 
layer  equation  is  integrated  from  bow  to  stern.)  Effectively,  calculations  along  a  series  of  streamlines  in 
the  three-dimensional  procedure  are  reduced  to  a  single  "average"  streamline,  namely  the  generating 
radius  R  =  R{x)  of  the  axisymmetric  body. 


R 


1 

O.FP 


Fig.  2  —  Axisymmetric  body  configuration  and  coordinate  system 


The  three  terms  on  the  right  hand  side  of  Eq.  (22)  show  that  the  boundary  layer  grows  due  to 
skin  friction,  adverse  pressure  gradient  ( dUjds  <  0),  and  streamline  convergence  ( dR/ds  <  0). 
Granville  suggests  that  the  pressure  distribution  be  calculated  by  using  the  Ry  equivalent  body,  and 
that  the  boundary  layer  growth  be  calculated  by  using  the  Rw  equivalent  body.  This  would  tend  to 
approximate  the  overall  friction  drag  acting  on  the  actual  hull  surface.  However,  it  is  not  clear  that  the 
Rw  body  gives  a  better  approximation  to  the  streamline  convergence  term  than  the  Rv  body.  Equa¬ 
tions  (20)  and  (21)  show  that  in  the  case  of  a  ship  with  a  vertical  stern  contour  of  zero  width,  the  case 
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for  most  of  the  ships  considered  in  the  present  study,  Rv  and  Rw  take  on  the  following  values  at  the 
stern  end 

^  =  0  (24a) 

and 

Rw=  —  (24b) 

77 

where  H  is  the  draft  of  the  actual  ship.  Consideration  of  the  last  term  in  Eq.  (22)  and  Eq.  (24)  shows 
that  the  Rv  body  would  tend  to  have  a  significantly  larger  streamline  convergence  effect  than  the  Rw 
body.  In  fact,  it  will  be  shown  later  (in  Fig.  4)  that  the  values  of  dR/ds  for  the  Rw  body  tend  to  have 
small  magnitudes,  giving  rise  to  weak  streamline  convergence  effects.  If  the  streamlines  on  the  hull  are 
traced  according  to  strict  potential  flow  theory,  they  would  all  converge  to  a  single  stagnation  point  at 
the  aft  perpendicular.  However,  due  to  viscous  flow,  this  stagnation  point  does  not  exist.  Also,  if  one 
uses  a  numerical  method  where  the  ship  surface  is  represented  by  a  series  of  panels,  the  approximate 
nature  of  the  calculation  leads  to  the  traced  streamlines  at  the  stern  showing  definite  convergence,  but 
not  to  a  point  [31-33].  The  actual  average  streamline  convergence  behavior  would  tend  to  lie  between 
those  corresponding  to  the  Rv  and  Rw  bodies.  Thus,  to  bracket  the  actual  effect,  boundary  layer  cal¬ 
culations  were  made  using  both  equivalent  bodies. 

Taking  the  distribution  of  Ue  to  be  known  (this  calculation  will  be  described  below),  two  addi¬ 
tional  equations  are  required  to  solve  for  the  3  unknowns  0,  Cy,  and  H  in  Eq.  (22)  The  two  additional 
equations  were  taken  to  be  the  well-known  Ludwieg-Tillmann  friction  relation  [34]  and  the  entrainment 
equation  of  Head  [35]  and  Standen  [36]: 


Cf  =  0.246  x  10-°'678"  Rg0m 

(25) 

-j-  ( 6G )  +  9(\^~  +  4  -7t|  =  F(G), 

ds  [Ue  ds  R  ds 

(26a) 

G  =  1.535  (H  -  0.7)"2-715  +  3.3, 

(26b) 

F(G)  =  0.0306  (G  -  3.0)~°'653 

(26c) 

where  Re  =  U9/v  is  the  Reynolds  number  based  on  momentum  thickness  (26d)  and  v  is  the  fluid 
kinematic  viscosity.  Equations  (22)  and  (26)  may  be  viewed  as  two  equations  for  the  two  unknowns  9 
and  //,  with  Eq.  (25)  furnishing  an  expression  for  Cf  which  appears  in  Eq.  (22). 

A  more  complex  formulation  for  Cf  suggested  by  Granville  [37]  and  used  in  Ref.  33  was  also 
tried  for  several  cases.  It  was  found  that  in  nearly  all  cases,  there  was  little  difference  between  the  two 
formulations  in  the  calculated  values  for  the  displacement  thickness  8*,  the  boundary  layer  variable  of 
principal  interest  in  the  present  study. 

Choice  of  Initial  Conditions 

An  important  and  nontrival  part  in  boundary  layer  calculations  is  the  choice  of  initial  point  s  and 
starting  values  for  the  dependent  variables  9  and  H .  For  cases  where  the  calculations  are  made  for  a 
particular  body  for  which  experimental  data  are  available,  the  convenient  procedure  is  to  use  measured 
values  of  9  and  H  at  a  given  initial  station  s0  to  start  the  calculations.  This  procedure  is  not  possible 
for  the  present  generalized  study  where  a  wide  range  of  ship  cases  are  considered.  A  second  possibility 
is  to  use  the  straightforward  approach  outlined  by  Garcia  and  Zazurca  [38].  They  give  a  simple  formula 
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for  the  growth  of  9 2  of  the  laminar  boundary  layer  in  the  bow  region  as  well  as  various  formulations  for 
the  location  of  transition,  after  which  the  turbulent  boundary  layer  formulation  of  Eqs.  (22),  (25),  and 
(26)  may  be  used.  This  approach,  however,  has  the  following  limitations.  In  most  model  tests,  tur¬ 
bulence  stimulation  devices  are  placed  in  the  bow  region.  While  this  fixes  the  location  of  transition 
from  laminar  to  turbulent  flow,  the  resulting  characteristics  of  the  turbulent  boundary  layer  are  a  com¬ 
plex  function  of  the  location,  size,  and  shape  of  the  device.  For  those  cases  where  no  stimulation  de¬ 
vices  are  used,  the  pressure  gradient  in  the  bow  area  often  leads  to  a  complex  flow  region  of  laminar 
separation  with  turbulent  reattachment,  which  is  not  modeled  by  the  above  approach. 

Consequently,  the  following,  approximately  uniform  approach  was  used  to  make  the  parametric 
runs  for  the  various  ship  cases  described  in  the  following  chapter.  The  location  of  transition  is  fixed  at 
a  value  of  Sq  corresponding  to  X/ L  =  0.05,  a  typical  location  for  turbulence  stimulators  used  in  model 
tests.  The  initial  value  of  the  momentum  thickness  0O  is  set  equal  to  0.7  of  the  corresponding  value  for 
a  flat  plate  with  the  same  length  Reynolds  number,  while  the  initial  value  of  the  shape  factor  H0  is  set 
equal  to  the  flat  plate  value,  usually  in  the  1.4  to  1.5  range.  This  was  largely  based  on  the  results  of 
Refs.  33,  39,  and  40.  These  results  show  that  the  boundary  layer  thicknesses  9  and  8*  in  the  bow 
region  are  usually  somewhat  smaller  than  flat  plate  values,  due  to  the  favorable  pressure  gradient  and 
streamline  divergence  in  this  region,  while  H  is  typically  1.4.  Numerical  experimentation  showed  that 
in  the  case  of  the  displacement  thickness  8*,  different  choices  of  9o  and  Hq  led  to  differences  in  the 
bow  region  which  tended  to  propagate  unchanged  to  the  stern  region.  These  differences  were  quite 
small  and  essentially  negligible  compared  to  values  of  8*  in  the  stern  region  where  most  of  the  bound¬ 
ary  layer  growth  takes  place  due  to  the  adverse  pressure  gradient  and  streamline  convergence. 


Calculation  of  Pressure  Distribution  on  Body  Surface 


In  the  present  thin-boundary  layer  approach,  where  the  pressure  is  assumed  to  be  constant 
throughout  its  thickness,  the  edge  velocity  Ue  appearing  in  Eq.  (22)  is  simply  related  to  the  pressure  p 
on  the  body  by  Bernoulli’s  equation. 


ue  =  ujr^cp 


(27) 


where  Cp  =  ( p~Po )/ 


is  the  pressure  coefficient  and  po  is  the  ambient  pressure  far  from  the 


body. 


There  are  well-known  numerical  procedures,  such  as  the  Hess  and  Smith  method  [41],  for  calcu¬ 
lating  the  potential  flow  pressure  distribution  on  arbitrary  bodies  of  revolution  to  any  desired  degree  of 
accuracy.  While  the  method  in  Ref.  41  is  straightforward,  a  number  of  numerical  steps  are  involved, 
such  as  matrix  inversion  and  calculation  of  elliptic  integrals  due  to  the  use  of  ring  sources  on  the  body 
surface.  In  the  present  study,  the  simple  and  direct  slender-body  theory,  as  outlined  by  Karamcheti 
[42],  is  used  to  calculate  the  pressure  distribution.  The  accuracy  of  the  method  is  improved  by  using 
two  modifications  which  are  described  later.  In  addition  to  its  obvious  simplicity,  the  present  choice  of 
calculation  methods  was  also  based  on  consideration  of  the  pressure  accuracy  requirements  at  both  ends 
of  the  body  (where  slender-body  theory  is  most  inaccurate).  The  approach  for  obtaining  initial 
conditions  for  the  boundary  layer  calculations  described  above  essentially  makes  it  unnecessary  to  have 
an  accurate  description  of  the  pressure  at  the  bow  end.  At  the  stern  end,  it  is  well-known  that  the 
potential  flow  pressure  distribution  is  substantially  modified  by  the  viscous  flow. 

In  conventional  slender-body  theory,  the  longitudinal  axis  of  the  body  is  divided  into  a  series  of 
N  -  1  line  segments  of  length  AXh  over  each  of  which  is  placed  a  source  of  uniform  density  q,,  where 
AXj  and  qt  may  vary  from  segment  to  segment.  The  potential  of  the  source  Gs  is  given  by 


r  V(I-^)2  +  R 2 
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where  £  is  a  source  point,  between  the  ends  of  the  z'th  line  segment.  The  strength  9,  is  directly 
by 


Qt 


U_ 

4tt 


dA 

dX 


UIttR 


Vi 


dRi 


dX 


given 

(29) 


'"V  I*  » 

where  A  =  ttR$  is  the  cross-sectional  area  of  the  body  (see  Eq.  (20)).  Upon  using  the  following  defi¬ 
nitions  for  the  axial  and  radial  velocities  on  the  body  surface  (X,R),  respectively  v*  and  vR,  due  to  a 
source  of  strength  q,  at  point  (£,,0), 


v*  = 


dX 


Qi 


and 


%/a-e,)2  +  i?2 

Qi 


(30a) 


ycr-^2, 

integrating  the  above  expressions  over  the  *th  source  segment  and  summing  over  the  N  —  1  source 
segments,  the  following  two  expressions  for  the  resultant  values  of  v*  and  vR  on  the  body  surface  due 
to  all  the  source  segments,  Vx  and  VR ,  are  obtained: 


(30b) 

source 


and 


vxa,R)  = 

U  N~l 

dA 

1 

1  1 

4*  k 

dX 

i 

J(X-£Li)2  +  R2 

y/(X-tni)2  +  R\ 

Vr(X,R)  = 

U 

dA 

X  ~  Zr, 

x  -  tu  1 

4nR  fy 

dX 

j  i 

ycr  -tRI)2  +  R2 

J(x-Ui)2  +  R2\ 

(31a) 


(31b) 


b  Ll  b  AM  - - -  w 

Noting  that  Vx/ U  is  a  higher  order  quantity  [42],  the  following  consistent  approximation  for  Ue 
on  the  body  surface  is  obtained  as 

t-  VRPW- 


Vx  Vi  vl 

1  +  2-V  +  1a+T W 


Vx  1  Vr 
1  H - -  +  ~ 


('37’) 


To  correct  the  inaccuracies  of  the  above  slender-body  theory  at  the  ends  of  the  body,  a  singularity 
gap  correction  proposed  by  Moran  [43]  is  applied  at  the  bow  and  stern  ends.  In  this  approach,  the 
source  distribution  does  not  extend  all  the  way  to  the  end,  but  there  is  a  gap  equal  in  length  to  one  half 
of  the  nose  radius  of  the  body.  The  gaps  at  the  bow  and  stern  ends,  respectively  yB  and  ys ,  then  are 
approximated  as 


Jb  = 


R NO 

2 


l 

4 


^2]  RHXJ  -  [fl2(Xi  =  0)  -  0] 

~dX  )*_0~  4X2 


(33a) 


and 


4 


dR2 

dX 


X-L 


RHXn  -  L)  -  RHXn.  1) 
4(Xn  -  XN-0 


(33b) 


where  R^o  and  Rnl  are  respectively  the  nose  radii  at  the  bow  and  stern  and  R  is  the  body  radius. 
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A  second  modification  is  made  to  account  for  the  effect  of  viscous  flow  on  the  stern  pressure  dis¬ 
tribution.  It  is  well-known  that  Ue  does  not  approach  zero  at  the  stern  end,  as  predicted  by  potential 
flow  theory.  The  modification  consists  of  an  iterative  procedure  similar  to  that  previously  used  by  the 
author  [44],  For  the  initial  boundary  layer  calculation,  the  calculated  values  of  Ue  over  the  last  5%  of 
the  body  are  discarded  and  are  replaced  by  a  linear  extrapolation  of  the  calculated  values  for  X/L  = 
0.95.  The  calculated  displacement  thickness  8*  is  then  added  to  the  body,  resulting  in  a  somewhat 
larger  body  with  a  stern  region  which  is  significantly  less  blunt.  For  the  second  (and  later)  boundary 
layer  calculations,  the  calculated  Ue  over  the  entire  displacement  body  (up  to  the  stern  end,  X/L  = 
1.0)  is  used.  This  is  possible  because  of  the  much  smoother  variation  of  the  shape  of  the  stern  region 
of  the  displacement  body.  The  above  procedure  is  repeated  for  a  total  of  four  iterations.  In  most 
cases,  there  is  little  difference  in  the  calculated  boundary  layer  characteristics  between  the  third  and 
fourth  iterations,  similar  to  the  trend  observed  in  Ref.  44.  However,  it  is  shown  later  that  the  method 
does  not  converge  satisfactorily  for  a  particular  case  of  a  ship  with  a  blunt  stern,  thus  giving  an  indica¬ 
tion  of  the  limit  of  applicability  of  the  present  procedure. 

Computerized  Calculation  Procedure 

The  above  calculations  have  been  implemented  in  the  form  of  a  computer  program  called  KBLST. 
The  program  first  accepts  input  for  the  half-width  /  of  the  ship  at  a  series  of  longitudinal  x  and  vertical 
y  stations.  The  input  values  of  /  may  be  from  a  data  file,  or  (as  is  the  case  for  most  of  the  ships  con¬ 
sidered  in  this  study)  they  may  be  analytically  defined.  The  program  then  calculates  the  wetted  area, 
the  normal,  and  the  source  strength  cr,  given  by  Eqs.  (15)  and  (17),  for  each  ship  panel  defined  by  the 
input  points.  At  each  input  x-station,  the  radii  Ry  and  R\y  in  Eqs.  (20)  and  (21)  for  the  equal  volume 
and  equal  wetted  area  bodies,  respectively,  are  calculated.  The  edge  velocity  equations  (Eqs.  (27)  to 
(33))  are  calculated  using  the  Rv  body,  while  the  boundary  layer  equations  (Eqs.  (22),  (23),  (25),  and 
(26))  are  calculated  using  either  the  Rv  or  Rw  body,  at  the  option  of  the  user.  If  viscous  flow  effects 
are  desired  in  the  wave  resistance  calculations,  the  final  calculated  displacement  thickness  is  added  to 
the  ship  hull  and  revised  values  of  the  source  strengths  <x  in  Eqs.  (15)  and  (17)  are  calculated.  In 
many  cases,  the  use  of  the  Rv  body  to  make  the  boundary  layer  calculations  gives  such  a  large  displace¬ 
ment  thickness  that  its  addition  to  the  original  hull  gives  a  resultant  displacement  hull  which  widens  as 
the  stern  end  is  approached,  i.e.,  gives  source  strengths  in  the  stern  region  which  are  opposite  in  sign  to 
those  for  the  original  hull.  It  was  felt  that  this  represented  an  unrealistically  large  modification.  In 
these  cases,  a  is  simply  set  equal  to  zero. 

The  wave  resistance  calculations  are  made  using  either  slender-ship  theory,  Eqs.  (15)  and  (16),  or 
thin-ship  theory,  Eqs.  (17)  and  (18).  The  principal  part  of  these  calculations  is  to  obtain  the  amplitude 
spectrum  as  a  function  of  longitudinal  wave  number  A2(kx)  by  using  equations  based  on  Eq.  (16)  in 
the  slender-ship  case  and  Eq.  (18)  in  the  thin-ship  case.  The  wave  resistance  is  then  calculated  by 
using  an  equation  equivalent  to  Eq.  (8).  The  program  also  contains  the  capability  of  calculating  the 
wave  elevation  over  a  rectangular  grid  on  the  free  surface  by  essentially  taking  the  inverse  Fourier 
transform  of  A2(kx). 

In  summary,  the  program  contains  six  options  for  calculating  wave  resistance  for  a  given  ship- 
Froude  number  case: 

1.  slender-ship— no  displacement  thickness, 

2.  slender-ship— displacement  thickness  using  Rw  body, 

3.  slender-ship— displacement  thickness  using  Ry  body, 

4.  thin-ship— no  displacement  thickness, 

5.  thin-ship— displacement  thickness  using  Rw  body, 

6.  thin-ship— displacement  thickness  using  Ry  body. 

For  most  of  the  ship  cases  considered  in  the  present  study,  all  six  options  were  exercised. 
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Due  to  the  simplicity  of  the  present  calculation  procedure,  the  program  KBLST  can  be  run  with 
ease  on  a  minicomputer  such  as  the  Hewlett-Packard  1000.  Running  time  for  the  usual  case  of  a  ship 
with  port  and  starboard  symmetry,  with  one  side  modeled  by  approximately  200  to  300  panels,  is  typi¬ 
cally  10  min  for  one  of  the  above  options.  This  makes  it  convenient  to  make  parametric  runs  for  the 
large  number  of  ship  forms  described  in  the  following  chapter. 

DESCRIPTION  OF  SHIP  HULLS 

Choice  of  Ship  Hulls 

The  choice  of  ship  hulls  was  based  on  a  number  of  considerations.  These  included  the  availability 
of  experimental  data  to  validate  the  theoretical  calculations,  the  existence  of  other  thin-  and  slender- 
ship  calculations  to  compare  with  the  present  approach,  the  need  to  cover  a  range  of  geometrical 
characteristics  to  illustrate  the  boundary  layer  effect,  and  the  existence  of  analytic  descriptions  for  hull 
shape  which  lead  to  convenient  and  precise  input  data. 

Based  on  these  considerations,  four  hulls  were  initially  chosen,  and  a  fifth  was  chosen  after  the 
start  of  the  calculations:  the  Wigley  hull  [5],  the  Series  60  (CB  =  0.60)  hull  [5],  the  Sharma  strut  [45], 
the  elliptic  bow-parabolic  stern  (EP)  strut  (B/ L  —  0.15,  where  B  —  ship  beam)  suggested  by  Noblesse 
[46],  and  a  thin  EP  (TEP)  strut  (B/L  =  0.05).  The  Wigley  and  Series  60  hulls  were  obvious  choices 
because  extensive  experimental  and  calculated  data  exist  for  these  hulls.  Both  are  among  the  test  hulls 
suggested  for  the  First  and  Second  DTNSRDC  (David  Taylor  Naval  Ship  Research  and  Development 
Center)  Workshops  on  Ship  Wave-Resistance  Computations.  The  Sharma  strut  represents  an  extremely 
thin  hull  form  for  which  a  series  of  experimental  data  and  calculations  using  thin-ship  theory  are  given 
in  Ref.  45.  No  experimental  data  exist  for  the  EP  and  TEP  struts.  However,  the  EP  strut  is  one  of  the 
suggested  hull  forms  for  calculation  for  the  Second  DTNSRDC  Workshop  and,  as  a  result,  a  series  of 
slender-ship  calculations  by  Scragg  [47]  are  available  for  comparison. 

Due  to  the  unsymmetrical  fore  and  aft  shapes  of  the  EP  strut,  it  was  originally  intended  to  make 
boundary  calculations  for  this  ship  moving  both  forward  and  backward.  These  would  show  if  the 
present  calculation  procedure  could  reproduce  the  general  trend  of  the  experimental  results  by  Wigley 
[48]  for  a  series  of  unsymmetrical  ships  towed  both  forward  and  backward.  Unfortunately,  the  ends  of 
this  hull  are  too  blunt  for  the  present  boundary  layer  calculation  approach  which  is  essentially  applicable 
for  slender-bodies.  This  was  demonstrated  by  the  fact  that  the  boundary  layer  calculations  did  not  con¬ 
verge.  For  example,  there  were  usually  large  differences  in  the  displacement  thickness  calculated  for 
the  third  and  fourth  iterations. 

Because  of  the  above  limitations,  the  TEP  hull  was  added.  The  TEP  hull  is  similar  to  the  EP  hull, 
the  difference  is  that  the  width  of  the  TEP  is  everywhere  only  1/3  that  of  the  EP.  The  boundary  layer 
calculations  converged  for  this  new  shape,  and  they  were  performed  for  the  ship  moving  both  forward 
(elliptic  bow)  and  backward  (parabolic  bow) . 

Geometric  Characteristics  of  Ship  Hulls 

Table  1  is  a  summary  of  the  principal  geometric  characteristics  of  the  hulls:  length  L,  beam-to- 
length  ratio  B/L,  beam-to-stern  length  ratio  B/Ls,  draft-to-length  ratio  H/ L,  block  coefficient  CB ,  and 
wetted  area  coefficient  C$,  where 

L  is  the  length  at  the  design  waterline, 

B  is  the  beam,  the  maximum  width, 

Ls  is  the  stern  length,  defined  as  the  axial  length  from  the  stern  end  to  the  point 
where  maximum  width  is  initially  attained, 

H  is  the  draft, 
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Table  1  —  Summary  of  Geometric  Characteristics 
of  the  Five  Chosen  Hulls 


Ship  Hull 

L 

(ft) 

B/L 

B/Ls 

H/L 

CB 

Cs 

Wigley 

20 

0.0625 

0.125 

0.100 

0.44 

0.66 

Series  60 

20 

0.133 

0.266 

0.053 

0.60 

0.71 

Sharma 

6.56 

0.05 

0.100 

0.150 

0.66 

0.95 

EP 

20 

0.15 

0.600 

0.075 

0.86 

0.95 

TEP 

20 

0.05 

0.200 

0.075 

0.86 

0.96 

C,  — 


LBH 
Sw 


is  the  block  coefficient, 

is  the  wetted  area  coefficient, 


L(2H  +  B ) 


VD  is  the  volume  of  the  underwater  hull,  and 
SV  is  the  wetted  surface  area. 


(34a) 

(34b) 


Note  that  for  all  the  hulls  considered  in  the  present  study,  the  same  definition  for  Ls  applied  to  the 
bow  area  gives  the  bow  length  LB  equal  to  Ls,  i.e., 

B/Lb  =  B/Ls.  (34c) 


The  table  shows  that  the  length  L  used  in  the  calculations  for  the  Wigley,  Series  60,  EP,  and  TEP  hulls 
is  6.1  m  (20  ft).  This  is  representative  of  the  length  used  in  model  tests  of  the  Wigley  and  Series  60 
hulls,  and  the  unsymmetrical  models  towed  in  both  the  forward  and  reverse  directions  in  Ref.  48.  A 
length  of  2  m  (6.56  ft)  is  used  for  the  Sharma  strut,  corresponding  to  that  used  in  the  model  tests  [45]. 

Table  2  gives  the  analytical  formulas  for  the  half-width  /  of  the  Wigley,  Sharma,  EP,  and  TEP 
hull  forms.  Figure  3  shows  the  stem  and  stern  contours  as  well  as  the  half-width  section  shapes  at  a 
series  of  longitudinal  stations  for  the  Series  60  hull  form. 

Figure  3  shows  that  the  Series  60  hull  has  a  nonvertical  stern  contour.  (In  the  present  study,  the 
small  deviations  of  the  stem  contour  from  the  vertical  are  neglected.)  The  other  four  hulls,  which  are 
analytically  defined,  have  vertical  stem  and  stern  contours.  Table  2  and  Fig.  3  show  that  the  Sharma, 
EP,  and  TEP  struts  are  wall-sided  (no  variation  of  half-width  /  with  depth  z)  while  the  Wigley  and 
Series  60  hulls  have  values  of  /  which  decrease  with  depth.  This  accounts  for  the  low  values  of  CB 
and  Cs  for  these  two  hulls.  The  Series  60  and  EP  hulls  are  significantly  wider  (larger  values  of  B/L) 
than  the  other  3  hulls.  The  stern  of  the  EP  hull  has  a  much  larger  value  of  B/Ls  than  the  other  four 
hulls  and  led  to  the  previously  mentioned  difficulty  in  the  boundary  layer  calculations.  This  in  turn  led 
to  the  necessity  of  adding  the  finer  TEP  hull,  which  has  a  stern  whose  value  of  B/ Ls  is  comparable  to 
those  of  the  remaining  hulls.  The  values  of  H/ L  range  from  the  shallow  draft  value  of  0.053  for  the 
Wigley  hull  to  0.15  for  the  deeply  submerged  Sharma  strut. 

The  Wigley  and  Sharma  hulls  are  symmetrical  fore  and  aft  while  the  remaining  three  hulls  are 
unsymmetrical.  The  EP  and  TEP  hulls  have  a  longitudinally  parallel  middle  body  section  of  constant 
width  amounting  to  one-half  of  ship  length.  The  remaining  three  hulls  have  widths  which  continually 
vary  with  x,  with  the  maximum  width  occurring  amidships. 
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7.6"  BOW  RADIUS 


SERIES  60  SECTION  CONTOURS 
Fig.  3  —  Geometric  description  of  Series  60,  CB  =  0.60  hull 
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Table  2  -  Formulas  for  the  Half- Width  / 
of  Four  Ship  Hulls 


*'=(*  +  1/2)/ (1/2) 

x  =  1  is  the  aft  perpendicular 

x'  =  +  1  is  the  forward  perpendicular 

1.  Wigley  hull 

f\  =  y  (1  -  *'2)[1  -  {z/H)2}  -  1  <  *'  <  +  1 

2.  Sharma  strut 


fi  =  y  (1  -  x'2) 


1  ^  x'  ^  +  1 


3.  Elliptic  bow-parabolic  (EP)  stern  strut 


f i  =  ~  [—  4x'(l  +  a:')] 

/3=  f 


f  -  B 


1  -  16 


2L 

2 


1/2 


"  1  *  *'  <  +  T 

+  y  ^  x'  <  +  1 


4.  Thin  elliptic  bow-parabolic  (TEP)  stern  strut 


/4=j/3 


~  1  <  X'  <  +1 


RESULTS 


Matrix  of  Computer  Runs 


Table  3  summarizes  information 
the  five  hulls.  The  table  gives: 


on  the  input  data  and  type  of  computer  runs  made  for  each  of 


NX,  the  number  of  stations  along  the  longitudinal  x-axis  at  which  the 
half-width  /  is  input  into  the  program, 

NZ,  the  number  of  stations  along  the  vertical  z-axis  at  which  /  is 
input, 

Fn  1 ,  the  lowest  Froude  number  Fn  for  which  wave  resistances  are 
calculated, 

Fn2,  the  highest  Froude  number  for  which  wave  resistances  are 
calculated, 

NF,  the  number  of  Fn  cases  for  which  wave  resistances  are  calculated, 

A/P,  the  number  of  calculations  for  each  Fn  case,  neglecting  boundary 
layer  effects,  and 

A/F,  the  number  of  calculations  for  each  Fn  case,  including  boundary 
layer  effects. 


Wioi/hQ  Va,Ue<n°f  NXa  tu  the  tlhree,hulls  which  have  continually  varying  widths  along  the  x-axis  (the 
igley,  Series  60,  and  Sharma  hulls)  is  25,  and  it  is  respectively  reduced  to  22  and  17  for  the  EP  and 
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Table  3  —  Summary  of  Computer  Runs 
for  the  Five  Ship  Hulls 


Ship  Hull 

NX 

NZ 

Fn  1 

Fn2 

NF 

NP 

NV 

Wigley 

25 

15 

0.16 

0.48 

11 

2 

4 

Series  60 

25 

7 

0.14 

0.36 

18 

2 

4 

Sharma 

25 

16 

0.21 

1.00 

15 

2 

4 

EP 

22 

16 

0.15 

1.00 

16 

2 

0 

TEP 

17 

16 

0.15 

1.00 

16 

2 

8 

TEP  hulls  which  have  parallel  middle  body  sections  extending  over  half  the  ship  length.  The  value  of 
NZ  is  15  or  16  for  the  four  hulls  which  are  analytically  defined  (the  Wigley,  Sharma,  EP,  and  TEP 
hulls).  It  is  reduced  to  7  for  the  Series  60,  which  corresponds  to  the  number  of  vertical  stations  at 
which  /  is  tabulated  in  Ref.  5.  The  Froude  number  range  for  the  three  hulls  for  which  model  test  data 
are  available  (the  Wigley,  Series  60,  and  Sharma  hulls)  correspond  to  the  experimental  range.  The 
Froude  number  range  for  the  EP  and  TEP  hulls  approximates  the  range  covered  by  the  calculations  of 
Scragg  [47]  for  the  EP  strut. 

Wave  resistance  calculations  neglecting  the  boundary  layer  effect  were  made  using  both  thin-ship 
and  slender-ship  theories  for  each  Froude  number  for  all  five  hulls.  For  the  Wigley,  Series  60,  and 
Sharma  hulls,  four  calculations  including  the  boundary  layer  displacement  thickness  8  *  were  also  made 
for  each  Froude  number:  thin-  and  slender-ship  theories,  each  using  8*  calculated  for  the  Ry  and  Rw 
bodies,  defined  in  Eqs.  (20)  and  (21),  respectively.  Due  to  the  previously  mentioned  nonconvergence 
of  the  boundary  layer  calculations,  no  wave  resistance  values  accounting  for  the  effect  of  8  *  are  shown 
for  the  EP  hull.  On  the  other  hand,  eight  such  calculations  were  made  for  each  Froude  number  for  the 
TEP  hull:  the  above-mentioned  four  calculations  for  the  hull  moving  both  forward  (blunter  elliptic 
bow,  finer  parabolic  stern)  and  backward  (finer  parabolic  bow,  blunter  elliptic  stern) . 


Calculated  Results 

Figure  4  shows  the  variation  of  Ry  and  Rwi  defined  respectively  in  Eqs.  (20)  and  (21),  along  the 
x-axis  for  the  equivalent  bodies  of  revolution  which  are  used  in  the  boundary  layer  calculations.  In  this 
figure,  the  and  Rv  curves  are  represented  by  solid  and  dotted  lines,  respectively.  Also,  the  Rv  and 
Rw  values  for  the  Sharma  strut  are  scaled  up  to  the  20-ft  ship  length  used  for  the  other  hulls,  recalling 
that  a  6.56-ft  length  is  actually  used  for  this  hull  (see  Table  1). 

Table  4  shows  the  ratio  8*/8*FP  for  the  Wigley,  Series  60,  Sharma,  and  TEP  hulls  for  five  Fn 
cases,  where  8*  is  the  displacement  thickness  for  the  Ry  and  R\y  bodies  at  stern  end,  and  8  PP  is  the 
corresponding  thickness  for  a  flat  plate  of  the  same  length  and  at  the  same  velocity.  For  the  TEP  hull, 
the  table  shows  this  ratio  for  the  hull  moving  both  forward  and  backward. 


Figure  5  shows  the  longitudinal  variation  of  the  ratio  (td/o-0  in  the  stern  region  for  the  hulls  con¬ 
sidered  in  Table  4  for  Fn  =  0.35,  where  o-0  and  aD  are  respectively  the  thin-ship  values  for  source 
strength  calculated  for  the  original  hull  and  the  hull  with  displacement  thickness  added.  As  mentioned 
previously,  if  the  displacement  thickness  modification  is  so  large  that  it  results  in  aD  having  a  sign 
opposite  to  that  of  o-0,  <td  is  simply  set  equal  to  zero.  The  present  approach  of  using  a  body  of  revolu¬ 
tion  results  in  a  constant  value  of  8  *  being  added  to  the  hull  at  each  longitudinal  station.  For  the  wall¬ 
sided  Sharma  and  TEP  struts,  this  results  at  a  given  x-station  in  a  value  of  crD/(To  which  does  not  vary 
with  depth  along  the  hull  and  is  the  value  shown  in  Fig.  5.  For  the  hulls  which  have  half-widths  which 
vary  with  depth,  the  Wigley  and  Series  60  hulls,  the  ratio  crD/cro  also  varies  with  depth.  For  these 
cases,  the  value  shown  in  the  figure  represents  an  average  over  depth. 
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Table  4  -  Ratio  8  */8  *FP  at  Stern  End  for  Equal- Wetted  Area 
(Rw)  and  Equal-Volume  (Rv)  Bodies  of  Revolution 


Wij 

s'ley 

Series  60 

Sharma  Strut 

Rw 

Rv 

Rw 

Ry 

Rw 

Ry 

Fn  =  0.2 

Fn  =  0.35 
Fn  =  0.5 

Fn  =  0.7 

Fn  =  1.0 

1.17 

1.15 

1.14 

1.13 

1.12 

5.29 

5.29 

5.27 

5.24 

5.22 

3.38 

3.79 

8.33 

8.44 

1.44 

1.42 

1.40 

1.38 

1.36 

5.64 

5.70 

5.74 

5.78 

5.81 

TEP  Strut 

Forward 

Backward 

Elliptic  Bow- 

Parabolic  Bow- 

Parabolic  Stern 

Elliptic  Stern 

Rw 

Ry 

Rw 

Ry 

Fn  =  0.2 

1.56 

5.04 

1.67 

3.34 

Fn  =  0.35 

1.54 

5.05 

1.66 

3.36 

Fn  --  0.5 

1.52 

5.05 

1.65 

3.37 

Fn  =  0.7 

1.50 

5.05 

1.64 

3.37 

Fn  =  1.0 

1.49 

5.04 

1.64 

3.38 

-20 


-19 
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Fig-  5  —  Variation  of  ratio  of  thin-ship  source  strengths  with  and 
without  viscous  modification  cr d/ctq  in  stern  region 
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Fieures  6  7  8  9,  and  10  show  the  calculated  values  of  the  wave  resistance  coefficient  Cw  as  a 
function  of  Fn  for  the  Wigley,  Series  60,  Sharma,  EP,  and  TEP  hulls,  where  Cw  follows  the  definition 
used  for  both  DTNSRDC  Workshops  on  Wave  Resistance: 


Cw  - 


D 


w 


\?U2SW 


(35) 


The  calculated  results  are  also  tabulated  in  Tables  5  to  9  for  the  above  hulls.  The  notation  for  the  cal¬ 
culated  results  in  the  figures  is  uniform  and  is  as  follows.  Points  plotted  by  the  symbols  , . x,  and  , 
denote,  respectively,  wave  resistance  values  with  no  8*  effect,  8  using  the  Rw  body,  and  8  us _  g 
Ry  body.  Values  obtained  by  using  slender-ship  theory  are  enclosed  by  a  circle,  while  those  obtai 
by  thin-ship  theory  are  left  uncircled.  Most  of  the  calculated  results  are  shown  in  these  figures.  How¬ 
ever  Aere  are  several  exceptions.  In  the  case  of  the  thin  Sharma  strut,  the  slender-  and  thin-ship 
values  (both  with  and  without  the  8*  effect)  are  so  close  that  the  slender-ship  results  have  been  omit¬ 
ted.  Similarly,  in  the  case  of  the  Series  60  hull  for  Fn  <  0.25,  the  thin-  and  slender-ship  results,  fo 
cases  where  the  8*  effect  is  included,  are  sufficiently  close  that  the  slender-ship  results  have  again  bee 
omitted.  In  these  cases,  the  reader  is  referred  to  the  corresponding  tables  for  a  more  detailed  look  at 

the  calculated  Cw  values. 

Fieures  6  7  and  8  show  experimental  data  for  the  Wigley,  Series  60,  and  Sharma  hulls.  A  range 
of  data  from  a  number  of  experiments  are  shown  for  the  Wigley  and  Series  60  hulls,  while  the  data 
from  the  single  series  of  experiments  reported  in  Ref.  45  are  shown  for  the  Sharma  strut.  The  majority 
of  the  experimental  data  for  the  Wigley  and  Series  60  hulls  reported  in  Ref.  5  are  for  the  case  where 
the  ship  model  is  free  to  sink  and  trim.  This  refers  to  the  fact  that  the  draft  and  pitch  angle  (rotation 
about  the  y  axis)  of  the  ship  hull  moving  at  forward  speed  are  generally  different  from  those  for _the 
ship  at  rest  These  changes  in  draft  (sinkage)  and  pitch  (trim)  arise  due  to  the  perturbation  flow  field 
induced  by  the  motion  of  the  ship.  The  numerical  results  of  Dawson,  Gadd,  and  others,  reported  in 
Ref.  5  indicate  that  the  effect  of  sinkage  and  trim  is  to  increase  the  wave  resistance  values  as  compared 
to  those  for  a  fixed  model.  Since  the  present  calculations  are  made  for  a  fixed  model,  it  is  desirable  to 
convert  the  experimental  data  to  this  case.  Fortunately,  Chen  and  Noblesse  [49]  have  performed  this 
conversion  for  the  Wigley  hull,  using  as  a  guide  the  reduction  factors  given  in  the  previously  men¬ 
tioned  numerical  results  [5],  Thus,  the  experimental  range  shown  in  Fig.  6  is  taken  from  Ref.  49  and 
applies  for  the  fixed  model  case.  Figure  6  also  shows,  as  discrete  points,  the  recent  experimental 
results  of  Kajitani  et  al.  [50]  specifically  for  a  fixed  model. 

The  experimental  range  shown  for  the  Series  60  in  Fig.  7  is  taken  from  Ref.  5,  and  the  results  are 
largely  for  cases  where  the  model  is  free  to  sink  and/or  trim.  Only  the  lowest  set  of  data  is  for  the 
fixed  model  case.  Finally,  the  experiments  for  the  Sharma  strut  [45]  are  for  the  fixed  model  case. 

Figures  6  to  9  also  show  calculations  of  other  authors  using  similar  thin-ship  and  slender-ship  cal¬ 
culations  to  compare  with  the  present  calculations.  Figure  6  shows  thin-  and  slender-ship  ca  cu  ations 
by  Lackenby,  on  p.  345  of  Ref.  5,  for  the  Wigley  hull.  Figures  7  and  9  show  slender-ship  calculations 
by  Scragg  [47]  for  the  Series  60  and  EP  hulls.  Figure  8  shows  thin-ship  calculations  [45]  for  the 
Sharma  strut.  In  all  of  these  calculations,  the  effect  of  viscosity  has  been  neglected. 


Figure  11  shows  the  differences  ACV  between  the  calculated  values  of  Cw 
forward  (finer  parabolic  stern)  and  backward  (blunter  elliptic  stern)  motions 

A  Cw  =  Cwp  ~  Cwb 


for  the  TEP  hull  in 
(36) 


where  CWF  and  CWB  are  respectively  the  values  of  Cw  for  forward  and  backward  motions, 
are  shown  for  thin-  and  slender-ship  theories  for  viscous  modifications  using  both  the 
bodies.  Table  10  shows  these  differences  in  tabulated  form. 


The  results 
R]y  and  Ry 
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Fig.  6  —  Wave  resistance  coefficients  for  Wigley  hull 
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Fig.  8  —  Wave  resistance  coefficients  for  Sharma  strut 
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Pig.  9  _  Wave  resistance  coefficients  for  EP  strut 
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Table  5  -  Calculated  Values  of  Cw  x  1000  for  Wigley  Hull 


Table  6  —  Calculated  Values  of  Cw  x  1000 
for  Series  60,  Cg  =  0.60  Hull 


8  * 

=  0 

8*  for  Rw  Body 

8  *  for  1 

Body 

Fn 

Thin 

Slender 

Thin 

Slender 

Thin 

Slender 

Ship 

Ship 

Ship 

Ship 

Ship 

Ship 

0.14 

1.221 

0.262 

0.114 

0.136 

0.162 

0.137 

0.15 

1.216 

0.304 

0.134 

0.163 

0.180 

0.161 

0.16 

1.100 

0.316 

0.173 

0.223 

0.273 

0.253 

0.17 

1.048 

0.286 

0.142 

0.199 

0.267 

0.257 

0.18 

1.135 

0.491 

0.295 

0.356 

0.373 

0.348 

0.19 

0.954 

0.479 

0.316 

0.402 

0.494 

0.465 

0.20 

1.267 

0.801 

0.560 

0.627 

0.634 

0.592 

0.21 

1.091 

0.620 

0.383 

0.446 

0.473 

0.430 

0.22 

1.017 

0.721 

0.555 

0.624 

0.748 

0.683 

0.23 

1.169 

0.926 

0.719 

0.805 

0.869 

0.828 

0.24 

1.046 

0.927 

0.657 

0.801 

0.837 

0.829 

0.25 

1.148 

1.217 

0.984 

1.113 

1.295 

1.203 

0.26 

1.987 

2.044 

1.955 

1.908 

2.308 

1.999 

0.28 

4.036 

3.398 

3.625 

3.095 

3.637 

2.954 

0.30 

3.898 

2.904 

3.162 

2.589 

2.932 

2.322 

0.32 

2.759 

2.163 

2.223 

1.951 

2.147 

1.789 

0.34 

2.779 

2.602 

2.578 

2.436 

2.713 

2.397 

0.36 

4.260 

4.020 

4.203 

3.847 

4.417 

3.841 
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Table  7  —  Calculated  Values  of  Cw  x  1000  for  Sharma  Strut 


Fn 

8* 

=  0 

8*  for  Rw  Body 

8  *  for  R  y  Body 

Thin 

Ship 

Slender 

Ship 

Thin 

Ship 

Slender 

Ship 

Thin 

Ship 

Slender 

Ship 

0.208 

0.208 

0.204 

0.145 

0.142 

0.092 

0.090 

0.219 

0.137 

0.135 

0.127 

0.124 

0.170 

0.166 

0.229 

0.317 

0.307 

0.294 

0.288 

0.287 

0.280 

0.243 

0.405 

0.396 

0.320 

0.313 

0.213 

0.207 

0.258 

0.225 

0.221 

0.193 

0.189 

0.207 

0.201 

0.276 

0.524 

0.514 

0.519 

0.509 

0.532 

0.521 

0.301 

0.922 

0.900 

0.810 

0.791 

0.589 

0.573 

0.334 

0.501 

0.488 

0.428 

0.417 

0.335 

0.325 

0.365 

0.667 

0.659 

0.677 

0.667 

0.760 

0.747 

0.408 

1.953 

1.922 

1.920 

1.888 

1.857 

1.823 

0.447 

2.937 

2.873 

2.814 

2.753 

2.505 

2.448 

0.500 

3.382 

3.282 

3.193 

3.101 

2.700 

2.614 

0.578 

3.022 

2.890 

2.833 

2.713 

2.318 

2.216 

0.709 

2.038 

1.883 

1.910 

1.767 

1.541 

1.421 

1.000 

0.875 

0.705 

0.823 

0.666 

0.661 

0.532 

Table  8  —  Calculated  Values  of 
Cyy  x  1000  for  (Elliptic 
Bow-Parabolic  Stern)  EP  Strut 


Fn 

8* 

=  0 

Thin 

Ship 

Slender 

Ship 

0.15 

6.61 

1.42 

0.20 

11.20 

3.95 

0.22 

8.97 

3.29 

0.24 

22.30 

8.26 

0.26 

15.97 

4.41 

0.28 

15.58 

8.46 

0.30 

27.68 

11.93 

0.35 

23.16 

7.86 

0.40 

14.45 

11.29 

0.45 

16.47 

15.57 

0.50 

17.87 

19.75 

0.60 

15.41 

30.95 

0.70 

11.79 

44.64 

0.80 

9.12 

62.34 

0.90 

7.06 

86.62 

1.00 

5.57 

115.06 
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Table  9  —  Calculated  Values  of  Cw  x  1000  for  (Thin 
Elliptic  Bow-Parabolic  Stern)  TEP  Strut  in  Forward  (Parabolic 
Stern)  and  Backward  (Elliptic  Stern)  Motion 


8*  = 

=  0 

Forward  Motion 

Backward  Motion 

8  *(RW) 

8  *(Ry) 

8*CM 

8*C Rv) 

Fn 

TS 

SS 

TS 

SS 

TS 

SS 

TS 

SS 

TS 

SS 

0.15 

1.14 

0.70 

1.10 

0.66 

1.24 

0.74 

0.83 

0.55 

0.39 

0.29 

0.20 

1.96 

1.35 

1.91 

1.30 

2.11 

1.42 

1.55 

1.17 

0.94 

0.81 

0.22 

1.52 

0.95 

1.45 

0.89 

1.52 

0.95 

1.20 

0.84 

0.74 

0.62 

0.24 

3.77 

2.87 

3.57 

2.69 

3.11 

2.22 

3.26 

2.65 

2.36 

2.11 

0.26 

2.64 

1.80 

2.44 

1.62 

1.93 

1.17 

2.25 

1.67 

1.59 

1.34 

0.28 

2.64 

1.95 

2.55 

1.86 

2.55 

1.79 

2.35 

1.87 

1.91 

1.68 

0.30 

4.63 

3.58 

4.41 

3.38 

3.81 

2.79 

4.19 

3.42 

3.40 

3.00 

0.35 

3.81 

2.59 

3.55 

2.39 

2.68 

1.65 

3.42 

2.50 

2.73 

2.20 

0.40 

2.39 

1.50 

2.26 

1.41 

1.90 

1.07 

2.18 

1.46 

1.82 

1.36 

0.45 

2.72 

1.74 

2.60 

1.66 

2.15 

1.27 

2.53 

1.68 

2.14 

1.59 

0.50 

2.94 

1.84 

2.82 

1.77 

2.23 

1.32 

2.75 

1.73 

2.33 

1.66 

0.60 

2.53 

1.52 

2.42 

1.49 

1.86 

1.14 

2.37 

1.31 

2.02 

1.26 

0.70 

1.93 

1.22 

1.85 

1.21 

1.41 

0.98 

1.82 

0.94 

1.56 

0.84 

0.80 

1.49 

1.04 

1.44 

1.06 

1.10 

0.93 

1.41 

0.72 

1.21 

0.56 

0.90 

1.16 

0.95 

1.11 

0.96 

0.84 

0.86 

1.09 

0.60 

0.94 

0.37 

1.00 

0.91 

0.75 

0.88 

0.80 

0.67 

0.81 

0.87 

0.41 

0.75 

0.20 

Table  10  —  Differences  in  Calculated  Values  of 
Cw  x  1000  for  TEP  Strut  in  Forward 
Motion  Minus  Corresponding  Value  in  Backward  Motion 


8*  for  Rw  Body 

8*  for  Rv  Body 

Fn 

Thin 

Ship 

Slender 

Ship 

Thin 

Ship 

Slender 

Ship 

0.15 

0.27 

0.11 

0.85 

0.45 

0.20 

0.36 

0.13 

1.17 

0.61 

0.22 

0.25 

0.05 

0.78 

0.33 

0.24 

0.31 

0.04 

0.75 

0.11 

0.26 

0.19 

-0.05 

0.34 

-0.17 

0.28 

0.20 

-0.01 

0.64 

0.11 

0.30 

0.22 

-0.04 

0.41 

-0.21 

0.35 

0.13 

-0.11 

-0.05 

-0.55 

0.40 

0.08 

-0.05 

0.08 

-0.29 

0.45 

0.07 

-0.02 

0.01 

-0.32 

0.50 

0.07 

0.04 

-0.10 

-0.34 

0.60 

0.05 

0.18 

-0.16 

-0.12 

0.70 

0.03 

0.27 

-0.15 

0.14 

0.80 

0.03 

0.34 

-0.11 

0.37 

0.90 

0.02 

0.36 

-0.10 

0.49 

1.00 

0.01 

0.39 

-0.08 

0.61 
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DISCUSSION  OF  RESULTS 

The  shapes  of  the  equivalent  bodies  of  revolution  are  discussed  first.  The  source  strength  modifi¬ 
cations  resulting  from  the  boundary  layer  calculations  using  these  bodies  are  discussed  next.  The  wave 
resistance  calculations,  in  terms  of  their  comparison  with  other  theoretical  values  as  well  as  with  experi¬ 
mental  data,  are  discussed  for  each  individual  hull.  Finally,  some  overall  trends  are  indicated. 

Equivalent  Bodies  of  Revolution 

Figure  4  shows  that  for  all  five  hulls,  the  equal-wetted  area  (R«,)  body  has  a  larger  radius  than 
the  corresponding  equal-volume  ( Rv )  body.  With  the  exception  of  the  stern  end  of  the  Series  60  hull, 
the  difference  between  the  two  bodies  is  largest  at  both  ends,  with  Rv  approaching  0  and 
\(l/Ry)dRyldX\  »  \(\lRw)dRwldX\.  This  leads  to  the  previously  noted,  larger  streamline  conver¬ 
gence  effect  (see  Eq.  (22))  for  the  Ry  body. 

The  exception  at  the  stern  end  in  the  case  of  the  Series  60  hull  is  simply  due  to  its  nonvertical 
stern  contour  (shown  in  Fig.  3)  effectively  resulting  in  the  ship  having  zero  draft  at  the  stern  end. 
According  to  Eq.  (24b),  in  this  case  is  also  0. 

The  asymmetry  of  the  EP  and  TEP  is  clearly  shown  by  the  fact  that  \dR/dx\  for  both  the  Ry  and 
Rw  bodies  is  larger  at  the  elliptic  bow  end. 

Boundary  Layer  Modification  to  Source  Strength 

Table  4  shows  that,  with  the  exception  of  the  Series  60  hull,  the  use  of  the  Rw  body  leads  to 
values  6*  which  are  only  1.1  to  1.7  times  the  corresponding  values  for  a  fiat  plate.  On  the  other  hand, 
the  use  of  the  Ry  body  gives  values  of  8*  which  are  3.4  to  8.4  times  the  flat  plate  values.  This  is 
largely  due  to  the  previously  mentioned  larger  streamline  convergence  effect  for  the  Ry  body.  In  the 
case  of  the  Series  60  hull,  which  has  a  nonvertical  stern  contour  resulting  in  Rw  =  0  at  the  stern  end, 
8*  for  the  Rw  body  is  3.8  times  the  flat  plate  value.  The  above  results  show  that  the  streamline  con¬ 
vergence  term  tends  to  dominate  the  growth  of  the  boundary  layer  in  the  stern  region.  Since  it  has  pre¬ 
viously  been  shown  that  the  Rw  and  Rv  bodies  tend  to  bracket  the  actual  streamline  convergence 
effect,  the  values  of  8*  calculated  on  these  two  bodies  may  also  tend  to  bracket  the  actual  value  of  8*. 
The  table  also  shows  that  over  a  fivefold  change  in  Froude  number,  which  (for  a  given  model  length) 
is  equivalent  to  a  fivefold  change  in  Reynolds  number,  there  is  little  change  in  the  ratio  8*/8*FP. 

Figure  5  shows  the  ratio  o-D/<r0  of  the  modified  and  original  values  of  source  strength  in  the  stern 
region  of  the  hulls.  In  accordance  with  the  previously  noted  fact  that  8*  for  the  Rw  body  is  always  less 
than  the  corresponding  value  for  the  Rv  body,  the  ratio  cr^/cro  using  8*  for  the  Rv  body,  at  a  given 
value  of  x,  is  always  less  than  that  for  the  R\y  case.  With  the  exception  of  the  case  for  the  TEP  hull 
moving  backward,  the  values  of  the  ratio  for  the  Ry  case  all  reach  zero  within  1  ft  of  the  stern  end.  On 
the  other  hand,  none  of  the  values  of  cr^/cro  reach  zero  for  the  Rw  case.  The  minimum  value  in  the 
case  of  the  Series  60  hull,  which  has  a  nonvertical  stern  contour  and  Rw  =  0  at  the  stern  end,  is 
approximately  0.1.  For  the  other  hulls,  the  minimum  value  ranges  between  0.55  and  0.90. 

Again,  with  the  exception  of  the  TEP  hull  moving  backward,  the  ratio  increases  with  increasing 
distance  from  the  stern  end.  At  x  =  -16,  which  corresponds  to  20%  of  L  away  from  the  stern  end,  all 
the  values  of  are  greater  than  0.9,  i.e.,  within  0.1  of  1.0,  which  is  the  value  corresponding  to  no 

modification  due  to  8*.  This  is  to  be  expected  because  8*  is  largest  at  the  stern  end  and  decreases 
away  from  the  stern. 

The  anomalies  in  the  case  of  the  TEP  hull  moving  backward  may  be  largely  explained  by  the 
bluntness  of  the  elliptic  end,  as  shown  in  Fig.  4,  and  the  thin-ship  definition  of  a  given  by  Eq.  (17). 
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Figure  4  shows  that  the  elliptic  end  of  TEP  (in  the  region  near  x  =  0)  is  clearly  the  second  bluntest 
end  (to  the  EP  elliptic  end  which  failed  to  converge  in  the  boundary  layer  calculations)  in  both  the  Ry 
and  Rw  cases.  The  figure  indicates  that  in  this  blunt  region  the  displacement  thickness  (which  is  added 
normal  to  the  hull  surface)  is  largely  in  the  x  direction,  with  a  relatively  small  contribution  to  the  thick¬ 
ness  /.  Equation  (17)  shows  that  the  resulting  value  of  a-D  would  be  relatively  close  to  <r0,  giving  the 
anomalous  increase  in  o-D/<r0  at  the  stern  end.  This  is  one  indication  of  the  inapplicability  of  thin-ship 
theory  for  blunt  ends. 

The  behavior  of  crD/cro  in  the  stern  region  may  be  compared  with  the  constant  reduction  factor  of 
0.6  for  the  entire  stern  used  by  Havelock  [181.  If  one  considers  the  region  nearest  the  stern  end 
(where  most  of  the  change  in  hull  geometry  takes  place)  as  most  significant  and  takes  the  point  5%  of 
L  away  from  the  stern  end  as  being  representative,  Fig.  5  shows  that  values  for  the  Ry  case  are 
clustered  around  0.4  while  those  for  the  Rw  case  range  from  0.65  to  0.9.  These  values  tend  to  bracket 
the  factor  assumed  in  Ref.  18. 


Wigley  Hull 

Figure  6  shows  values  of  wave  resistance  coefficients  calculated  by  the  present  approach,  the 
theoretical  results  by  Lackenby  as  reported  in  Ref.  5,  the  range  of  experimental  data  corrected  for  sink- 
age  and  trim  [49],  and  the  recent  measured  data  [50]  specifically  for  a  fixed  model. 

Perhaps  the  most  noticeable  feature  of  the  figure  is  the  striking  agreement  between  the  thin-ship 
and  slender-ship  results,  neglecting  viscous  effects,  of  the  present  approach  and  those  of  Lackenby. 
The  agreement  is  usually  to  within  plotting  accuracy.  This  agreement  confirms  the  correctness  of  the 
numerical  part  of  the  present  potential  flow  approaches. 

Without  viscous  flow  effects,  the  thin-ship  values  are  always  higher  than  corresponding  slender- 
ship  values.  Similarly,  for  a  given  boundary  layer  approach  using  either  the  Rw  or  Rv  body  the  thin- 
ship  values  are  again  higher  than  those  using  slender-ship  theory. 

As  in  previous  studies,  the  effect  of  the  boundary  layer  is  usually  to  reduce  the  value  of  Cw.  For 
both  slender-  and  thin-ship  theories,  the  reduction  is  larger  for  the  approach  using  the  Ry  body  than 
for  the  Rfy  body.  This  is  due  to  the  larger  values  of  8*  for  the  Ry  body.  For  a  given  ship  theory,  the 
difference  in  Cw  between  the  two  boundary  layer  approaches  ranges  from  nearly  zero  at  certain  Fn  to 
0.4  x  10~3  at  other  Fn.  The  reduction  for  the  approach  using  the  Rw  body  is  always  quite  small. 

For  Fn  <  0.2,  all  the  calculated  values  lie  above  the  experimental  range.  For  Fn  >  0.2,  consid¬ 
ering  the  Cw  values  for  both  the  Ry  and  Rw  bodies,  the  calculated  values  for  thin-ship  theory  tend  to 
either  lie  within  the  experimental  range  or  bracket  the  upper  limit.  The  calculated  values  for  slender- 
ship  theory  tend  to  either  lie  within  the  experimental  range  or  bracket  the  lower  limit.  The  recent 
experimental  data  specifically  for  a  fixed  model  [50]  consistently  fall  near  the  lower  limit  of  the  experi¬ 
mental  range  [49].  Noting  that  the  results  in  Ref.  49  are  obtained  by  applying  correction  factors  to  data 
originally  obtained  for  models  free  to  sink  and  trim,  the  data  in  Ref.  50  tend  to  indicate  that  the  lower 
range  of  the  experimental  data  in  Ref.  49  may  be  more  appropriate  for  a  fixed  model. 

On  the  whole,  the  trend  of  the  viscous  correction  is  encouraging.  If  the  boundary  layer  calcula¬ 
tions  for  the  Rw  and  Ry  bodies  are  viewed  as  bracketing  the  actual  viscous  flow  effect,  the  results  indi¬ 
cate  that  a  more  accurate  boundary  layer  modeling  will  bring  the  thin-ship  calculations  into  closer  agree¬ 
ment  with  experimental  data  and  the  slender-ship  calculations  into  rather  good  agreement  for  most  Fn 
cases. 


32 


NRL  REPORT 8881 


Series  60  Hull 
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visddus  efffecfS3,  fdr'the  present  apprdachare  in  reasonably 
good  bufnot  pretise,  a^reemertf  with’thbSe  5f  Ref.  47;1*  Quadrilateral  panels  are%sed  in  the  present 
approach  While  triangular  panfels'  ar^uSed3in  Ref.  47V3Tt  wiiDbe'shownlater  in  the  case  Of  the  EP  hull 
where  quadrilateral  panels  are  used  in  both  approaches  that  the  agreement  is  to  within  plotting  accuracy. 

As  in  the  case  of  the  Wigiey  hull  in  the  absence  of  viscous  effects,  the  thin-ship  values  are  usually 
higher  than  corresponding  slender-ship  values.  When  viscous  effects  . are  included,  some  trends  are  the 
same  as  those  for  the  Wigiey  hull  while  others  are  not.  Figure  7  and  Table  6  show  that  the  effect  of 
adding  8  *  to  the  hull  is  again  to  generally  (but  not  always)  reduce  the  value  of  Cw.  Unlike  the  Wigiey 
hull  case,  the  reduction  due  to  8*  calculated  for  the  Rw  body  may  be  large  and  often  slightly  exceeds 
the  reduction  for  the  Rv  case.  This  is  due  to  the  previously  mentioned  fact  that  Rw  also  goes  to  zero 
at  the  stern  end  for  this  hull,  which  has  a  nonvertical  stern  contour.  In  the  Rv  case,  the  thin-ship 
values  are  higher  than  the  slender-ship  values,  as  in  the  case  when  viscous  effects  are  neglected.  On 
the  other  hand,  in  the  Rw  case  the  slender-ship  values  are  higher  for  Fn  <  0.25. 


The  agreement  of  the  calculated  results  with  experimental  data  is  considerably  worse  than  for  the 
Wigiey  hull.  Even  with  the  boundary  layer  correction,  the  thin-ship  results  for  Fn  ^  0.25  generally  lie 
well  above  the  upper  limit  of  the  experimental  data.  In  this  Froude  number  range,  the  slender-ship 
results  (with  the  8  *  correction)  are  lower  than  the  thin-ship  results.  At  times  they  lie  within  the  exper¬ 
imental  range  for  the  free  model  case  but  are  considerably  higher  than  the  data  for  the  Fixed  model 
case.  For  the  lower  Fn ,  most  of  the  calculated  data,  with  the  notable  exception  of  the  potential  flow 
thin-ship  results,  lie  within  the  range  of  the  experimental  data  but  again  appear  to  be  higher  than  the 
extrapolation  of  the  fixed  model  data.  It  appears  that  for  this  full  hull  form  (B/L  =  0.13)  with  a  fiat 
bottom,  a  more  complex  ship  theory  is  required.  The  results  shown  in  Ref.  5  indicate  that  the  use  of 
low  speed  theory,  such  as  that  by  Dawson  [8],  or  Guilloton’s  method  [11,12]  may  give  substantially 

more  accurate  results. 


Sharma  Strut 

As  mentioned  previously,  the  slender-ship  results  are  not  plotted  in  Fig.  8  due  to^their  closeness 
to  the  thin-ship  results.  This  is  not  surprising  because  of  the  fineness  of  the  hull  (B/L  =  0.05).  Table 
7  shows  that  in  all  cases,  with  and  without  the  boundary  layer  modification,  the  values  of  Cw  calculated 
by  thin-ship  theory  are  slightly  higher  than  the  corresponding  slender-ship  values.  This  is  in  accordance 
with  the  trend  previously  noted  for  the  Wigiey  and  Series  60  hulls.  Also,  as  in  the  Wigiey  case,  the 
thin-ship  results  (with  no  8  *  effect)  calculated  by  the  present  approach  are  in  excellent  agreement  with 
corresponding  results  by  Sharma  [45]. 

As  in  the  case  of  the  previous  two  hulls,  the  effect  of  adding  8  *  is  to  reduce  the  value  of  In 
the  case  of  the  Wigiey  hull,  the  reduction  is  generally  small  for  the  Rw  case  and  of  variable  size 
(dependent  on  Froude  number)  for  the  Rv  case.  The  calculated  values  with  no  8  *  effect  are  already  in 
reasonably  good  agreement  with  experimental  results,  and  the  addition  of  the  small  reduction, for  the 
Rw  case  improves  the  agreement  at  the  high  Fn  range.  The  use  of  8*  for  the  Ry  case  generally 
worsens  the  agreement.  This  is  probably  because  for  this  deeply  submerged  hull  (H/L  =  0.15)  the 
use  of  two-dimensional  boundary  layer  theory,,  instead  of  the  present  axisymmetric  approach,  may  be 
more  appropriate  over  most  of  the  hull.  In  the  two-dimensional  approach,  where  the  streamlines  are 
assumed  not  to  change  in  the  z-direction,  the  streamline  convergence  term  in  Eq.  (22)  goes  to  zero. 
This  is  the  effect  which  is  emphasized  more  by  use  of  the  Ry  body  than  the  Rw  body. 
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EP  Strut  (Elliptic  Bow-Parabolic  Stern,  B/L  =  0.15) 

Figure  9  shows  the  thin-ship  and  slender-ship  values  calculated  by  the  present  approach  and  the 
slender-ship  calculations  using  a  similar  theory  by  Scragg  [47],  As  mentioned  previously,  the  boundary 
layer  calculations  for  this  hull  failed  to  converge  due  to  the  extremely  blunt  stern,  B/Ls  =  0.60. 

....  The  figure  shows  that’  excePt  for  the  cross-over  region  of  0.4  <  Fn  <  0.5,  there  are  usually  large 
differences  between  the  thin-ship  and  slender-ship  values.  This  is  not  surprising  because  of  the  fullness 
of  the  entire  hull,  B/L  =0.15  and  the  bluntness  of  the  bow  and  stern,  both  of  which  have  B/L  = 
0.6.  Unlike  the  other  hulls,  the  slender-ship  results  are  not  generally  lower  than  the  thin-ship  values 
They  are  substantially  higher  for  Fn  ^  0.5. 

The  slender  ship  values  of  Ref.  47  agree  to  within  plotting  accuracy  with  the  present  results.  This 
is  not  surprising  because  the  similarity  of  the  approaches  and  the  use  of  quadrilateral  panels  in  both 
studies.  This  contrasts  with  the  previous  Series  60  case,  where  the  use  of  triangular  panels  in  Ref.  47 
leads  to  small,  but  noticeable,  differences  between  the  two  approaches. 


TEP  Strut  (Thin  Elliptic  Bow-Parabolic  Stern,  B/L  =  0.05) 

Figure  10  shows  the  thin-  and  slender-ship  values  calculated  by  the  present  approach,  with  and 
without  viscous  effects,  for  the  ship  moving  both  forward  (finer  parabolic  stern)  and  backward  (blunter 
elliptic  stern).  Figure  11  shows  the  difference  A Cw,  Cw  for  forward  motion  minus  Cw  for  backward 
motion,  for  thin-  and  slender-ship  theories  using  5*  calculated  for  the  Rv  and  Rw  bodies. 

Note  that,  without  viscous  effects,  the  present  thin-  and  slender-ship  theories  should  give  the 
same  value  of  Cw  for  both  forward  and  backward  motions.  This  was  confirmed  by  several  computer 
runs.  Noblesse  [51]  points  out  that  this  need  not  be  the  case  in  higher  order  slender-ship  theories. 

Many  of  the  trends  shown  in  Fig.  10  are  similar  to  those  for  the  Wigley  hull.  Without  viscous 
effects,  the  thin-ship  values  are  always  higher  than  the  corresponding  slender-ship  values.  Similarly, 
for  a  given  direction  of  motion  and  a  given  8*  approach,  the  thin-ship  values  are  always  higher.  As 
before,  the  reduction  in  Cw  for  the  Rw  case  using  thin-ship  theory  is  small  for  forward  motion.  How¬ 
ever,  the  reduction  is  somewhat  larger  for  backward  motion  (blunt  stern)  in  certain  Fn  cases,  for  both 
thin-  and  slender-ship  theories.  This  is  because  8  */8  *FP  for  the  Rw  case  for  the  Wigley  hull  is  only  1  2 
and  increases  to  1.6-1. 7  for  the  TEP  hull. 

As  mentioned  previously,  one  of  the  major  points  of  interest  for  this  unsymmetrical  hull  is  to  see 
if  the  present  calculation  procedure  can  reproduce  the  experimental  trend  for  A  Cw  obtained  by  Wigley 
[48]  for  several  unsymmetrical  hulls  towed  in  both  directions.  The  general  experimental  trend  is  for 
Cw  to  be  larger  for  the  finer  stern,  i.e.,  ACV  >  0  in  the  range  0.2  <  Fn  <  0.5,  with  t±Cw  reaching  a 
maximum  at  Fn  —  0.3.  For  Fn  >  0.5,  the  experimental  results  show  little  difference  in  the  value  of 
Cw  in  either  direction  of  motion,  i.e.,  ACV  ~  0. 

Figure  1 1  shows  that  this  trend  is  reproduced  quite  well  by  the  thin-ship  approach  using  8  *  calcu¬ 
lated  for  the  Rw  body  and  approximately  for  the  Rv  body.  Unfortunately,  the  experimental  trend  is 
not  reproduced  by  the  slender-ship  results.  For  a  given  approach,  the  slender-ship  curve  is  approxi¬ 
mately  parallel  to,  but  below,  the  corresponding  thin-ship  curve  for  Fn  <  0.40.  This  often  results  in 
Cw  having  the  incorrect  sign.  For  Fn  >  0.5,  the  slender-ship  curve  increases  to  large  positive  values 
instead  of  approaching  the  nearly  zero  values  obtained  experimentally.  This  seems  to  indicate  that  the 
inaccuracies  in  the  present  boundary  layer  approach  are  magnified  in  the  case  of  slender-ship  theory.  In 
this  connection,  it  is  of  interest  to  note  that  in  most  previous  studies,  the  boundary  layer  modifications 
have  been  used  in  conjunction  with  thin-ship  theory. 
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Overall  Trends 

Except  for  the  EP  strut  (B/L  =  0.15),  which  is  the  fullest  ship  hull  considered  in  the  present 
study,  the  thin-ship  values  of  Cw  are  generally  higher  than  the  corresponding  slender-ship  values,  both 
with  and  without  boundary  layer  effects.  In  some  cases,  the  differences  are  quite  small  as  in  all  the  cal¬ 
culations  for  the  thin  Sharma  strut  (B/L  =  0.05)  or  at  the  low  Fn  range  for  the  Series  60  hull  results 
including  boundary  layer  effects.  In  the  case  of  the  EP  strut,  the  slender-ship  values  are  considerably 
higher  than  corresponding  thin-ship  values  for  Fn  >  0.5. 

The  addition  of  8  *  serves  generally,  but  not  always,  to  reduce  the  value  of  Cw.  Due  to  the  larger 
streamline  convergence  effect  for  the  Rv  body,  the  reduction  in  Cw  is  generally  larger  for  8*  calculated 
using  this  body  than  for  the  Rw  body.  The  exception  occurs  for  the  Series  60  hull,  which  has  a  non¬ 
vertical  stern  contour;  this  results  in  Rw  =  0  at  the  stern  end. 

The  values  calculated  by  the  present  approach,  without  viscous  effects,  compare  quite  well  with 
other  similar  approaches.  The  minor  discrepancy  in  the  case  of  the  Series  60  hull  is  due  to  the  use  of 
quadrilateral  panels  in  the  present  study  and  triangular  panels  in  Ref.  47.  Also,  as  required  by  theoreti¬ 
cal  considerations,  the  same  value  of  Cw  is  obtained  for  the  unsymmetrical  TEP  hull  moving  forward 
and  backward  in  the  absence  of  viscous  effects. 

Agreement  with  experimental  data  is  erratic  and,  as  may  be  expected,  is  largely  dependent  on  hull 
shape.  Generally,  the  slender-ship  results  are  in  relatively  better  agreement  with  experimental  data 
than  thin-ship  values.  Also,  the  addition  of  the  boundary  layer  displacement  thickness  generally 
improves  the  agreement.  In  the  case  of  the  thin  Sharma  strut  (B/L  =  0.05),  the  theoretical  calcula¬ 
tions  without  viscous  effects  are  in  reasonably  good  agreement  with  experimental  data.  The  addition  of 
8*  for  the  Rw  body  tends  to  improve  the  agreement,  while  8*  for  the  Rv  body  tends  to  give  too  large 
a  reduction  in  Cw  and  thus  worsens  the  agreement.  This  may  be  due  to  the  fact  that  the  use  of  the  Rv 
body  tends  to  exaggerate  the  streamline  convergence  effect  which  may  be  largely  absent  on  this  deeply 
submerged  strut.  In  the  case  of  the  Wigley  hull,  considering  8*  for  both  the  Ryy  and  Rv  cases,  the 
thin-ship  values  tend  to  either  lie  within  the  experimental  range  or  bracket  the  upper  limit;  the 
corresponding  slender-ship  values  either  lie  within  the  range  or  bracket  the  lower  limit,  which  according 
to  the  data  of  Ref.  50,  may  be  more  applicable  for  the  present  case  of  a  fixed  model. 

In  the  case  of  the  full  Series  60  hull  (B/L  =  0.133),  which  also  has  a  flat  bottom,  the  calculated 
results  generally  are  much  higher  than  the  measured  data  for  a  fixed  model.  For  this  hull,  it  appears 
that  the  boundary  layer  correction  alone  is  not  sufficient  to  obtain  reasonable  agreement.  Instead,  the 
use  of  more  complex  wave  resistance  theories,  such  as  low-speed  theory  or  Guilloton’s  method,  give 
improved  agreement. 

In  the  case  of  the  unsymmetrical  TEP  hull,  the  calculated  thin-ship  results  for  both  the  Rw  and 
Rv  cases  generally  reproduce  the  experimental  trend  (obtained  for  other  unsymmetrical  hulls)  for  the 
difference  in  Cw  between  forward  and  backward  motions.  On  the  other  hand,  the  slender-ship  results 
do  not  reproduce  this  trend.  This  seems  to  indicate  that  inaccuracies  in  the  boundary  layer  modeling 
tend  to  be  magnified  in  the  slender-ship  case. 

SUMMARY  AND  CONCLUSION 

Calculation  Approach 

A  relatively  simple  calculation  procedure  has  been  developed  for  obtaining  the  wave  resistance, 
including  viscous  effects,  on  monohull  ship  forms.  The  procedure  has  been  implemented  as  a  com¬ 
puter  program  which  runs  well  on  a  minicomputer,  such  as  the  Hewlett-Packard  1000. 
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The  wave  resistance  calculations  are  made  by  using  the  direct  thin-ship  and  zeroth  order  slender- 
ship  theories.  In  these  approaches,  the  source  strengths  are  obtained  directly  from  input  hull  geometry 
and  do  not  require  iterative  procedures.  The  modification  for  viscous  flow  is  obtained  by  adding  the 
boundary  layer  displacement  thickness  8  *  to  the  actual  ship  hull.  8  *  is  calculated  by  approximating  a 
double  model  of  the  actual  hull  by  two  equivalent  axisymmetric  bodies  which  respectively  have  volume 
(Rv  body)  and  wetted  area  ( Rw  body)  equal  to  those  of  the  double  model.  Such  an  axisymmetric 
approach  for  making  the  boundary  layer  modification  does  not  appear  to  have  been  previously  tried  in 
studies  which  model  the  effect  of  viscosity  on  wave  resistance.  The  standard  momentum  integral  equa¬ 
tion  is  used  to  model  the  boundary  layer  growth  along  the  axisymmetric  bodies.  It  appears  that  8  *  cal¬ 
culated  for  the  Rv  and  Rw  bodies  may  bracket  the  actual  boundary  layer  effect  since  they  respectively 
tend  to  overemphasize  and  underemphasize  the  streamline  convergence  term  which  generally  dom¬ 
inates  boundary  layer  growth  in  the  stern  region.  The  pressure  gradient  on  the  body  is  calculated  by 
using  ordinary  slender-body  theory,  modified  by  singularity  gaps  at  the  ends  dependent  on  body  nose 
radius  and  a  simple  iterative  procedure  to  account  for  the  boundary  layer  thickness  in  the  stern  region. 


Hull  Selection 

The  above  procedure  was  used  to  calculate  the  wave  resistance  coefficient  Cw  over  a  range  of 
Froude  numbers  Fn  for  five  ship  hulls:  Wigley,  Series  60  (CB  =  0.60),  Sharma,  (elliptic  bow-parabolic 
stern)  EP,  and  (thin  EP)  TEP.  These  hulls  range  in  fullness  ratio  B/L  from  0.05  to  0.15  and  draft 
ratio  H/L  from  0.053  to  0.15.  The  first  three  hulls  have  extensive  experimental  data.  The  fourth  was 
proposed  for  calculation  for  the  Second  DTNSRDC  Workshop  on  Wave  Resistance.  The  fifth  hull  was 
used  mainly  to  indirectly  compare  the  calculated  trend  in  the  difference  between  Cw  for  the  hull  in  for¬ 
ward  and  backward  motions  with  experimental  data  for  other  unsymmetrical  hulls. 

Six  calculations  of  Cw  were  made  at  each  Fn  for  the  Wigley,  Series  60,  and  Sharma  hulls:  thin- 
and  slender-ship  theories,  each  for  the  cases  of  no  8*  modification  and  8*  calculated  for  the  Rw  and 
Ry  bodies.  For  the  EP  hull,  only  results  with  no  viscous  effects  are  shown  since  the  boundary  layer 
calculations  failed  to  converge  for  this  form  which  is  the  fullest  of  the  five.  Ten  calculations  of  Cw 
were  made  at  each  Fn  for  the  TEP  strut:  thin-  and  slender-ship  theories  for  the  case  of  no  viscous 
effects,  and  these  two  theories  for  8*  using  the  Rw  and  Rv  bodies  for  motion  in  the  forward  as  well  as 
backward  directions. 

Calculated  Results 

With  the  exception  of  the  Series  60  hull  which  has  a  nonvertical  stern  contour,  the  value  of  8  *  at 
the  stern  end  of  the  Rw  axisymmetric  body  ranges  from  1.1  to  1.7  times  the  flat  plate  value  b*FP.  For 
the  Rv  axisymmetric  body,  the  calculated  value  of  8*  is  substantially  larger,  ranging  from  3.4  to  8.4 
times  8  }P.  In  many  cases,  the  addition  of  this  large  value  of  8*  results  in  a  hull  form  which  widens  as 
the  stern  end  is  reached,  which  would  give  a  source  strength  which  has  a  sign  opposite  to  that  of  origi¬ 
nal  hull.  In  these  cases,  the  new  source  strength  is  simply  set  equal  to  zero. 

Many  of  the  calculated  trends  are  similar  to  those  of  previous  studies.  For  example,  in  the 
absence  of  viscous  effects,  the  thin-ship  values  of  Cw  are  usually  higher  than  corresponding  slender- 
ship  values.  The  notable  exception  occurs  in  the  case  of  the  fullest  form  EP  hull,  for  which  the 
slender-ship  values  are  substantially  higher  for  Fn  >  0.5.  This  trend  also  generally  holds  when  the  8  * 
effect  is  included;  the  effect  of  adding  8*  is  to  generally,  but  not  always,  reduce  the  value  of  Cw.  With 
the  principal  exception  of  the  previously  mentioned  case  of  the  Series  60  hull,  the  reduction  in  the  Rw 
case  is  usually  smaller  than  that  in  the  Rv  case. 

In  the  absence  of  viscous  effects,  the  values  calculated  by  the  present  approach  are  in  excellent 
agreement  with  values  calculated  in  other  studies  using  comparable  theories,  for  the  Wigley,  Sharma, 
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and  EP  hulls.  The  minor  discrepancy  in  the  case  of  the  Series  60  hull  is  largely  due  to  the  use  of  quad¬ 
rilateral  panels  in  the  present  study  and  the  use  of  triangular  panels  in  the  other  study.  Also  in  the 
absence  of  viscous  effects,  the  present  approach  predicts  the  same  values  of  Cw  for  the  unsymmetrical 
TEP  strut  in  forward  and  backward  motions,  as  theoretically  required. 

In  terms  of  agreement  with  experimental  results,  the  slender-ship  values  are  generally  better  than 
corresponding  thin-ship  values.  Also,  for  either  ship  theory,  the  addition  of  8  *  generally  improves  the 
agreement.  The  agreement  with  experimental  data  for  the  individual  hulls  is  erratic  and,  as  may  be 
expected,  is  dependent  on  hull  shape.  The  agreement  is  generally  good  for  the  thin  Sharma  strut.  The 
use  of  6*  for  the  Rv  body,  which  emphasizes  the  streamline  convergence  effect,  somewhat  worsens  the 
agreement  since  this  effect  is  relatively  small  for  this  deeply  submerged  hull  for  which  a  two- 
dimensional  boundary  layer  approach  may  be  more  appropriate.  For  the  Wigley  hull,  the  use  of 
slender-ship  theory  with  8*  for  the  Rw  and  Rv  bodies  gives  results  which  tend  to  bracket  the  experi¬ 
mental  values.  Thus,  for  this  hull,  the  use  of  a  more  accurate  boundary  layer  model  may  give  good 
agreement  with  experimental  data  for  most  Fn  cases.  The  agreement  in  the  case  of  the  Series  60  hull  is 
quite  poor.  Even  in  the  case  of  the  slender  ship  with  8  *  modification,  the  calculated  values  are  usually 
considerably  higher  than  the  appropriate  experimental  data.  Results  given  at  the  First  DTNSRDC 
Workshop  on  Wave  Resistance  indicate  that  the  use  of  a  more  complex  wave  resistance  theory,  such  as 
low-speed  theory  or  Guilloton’s  method,  brings  substantially  improved  agreement. 

In  the  case  of  the  TEP  strut,  the  thin-ship  results  using  the  8*  correction  for  either  the  Rv  or  Rw 
body  generally  reproduce  the  experimentally  observed  trend  (for  other  unsymmetrical  hulls)  for  the 
difference  between  Cw  for  forward  and  backward  motions.  This  trend  is  not  reproduced  by  the 
slender-ship  calculations.  This  seems  to  indicate  that  slender-ship  theory  is  more  sensitive  to  inaccura¬ 
cies  in  the  boundary-layer  model  than  thin-ship  theory. 
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